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Abstract 



The identity of dark matter is one of the key outstanding problems in both 
particle and astrophysics. In this thesis, I describe a number of comple- 
mentary searches for particle dark matter. I discuss how the impact of 
dark matter on stars can constrain its interaction with nuclei, focussing 
on main sequence stars close to the Galactic Centre, and on the first stars 
as seen through the upcoming James Webb Space Telescope. The mass 
and annihilation cross-section of dark matter particles can be probed with 
searches for gamma rays produced in astronomical targets. Dwarf galax- 
ies and ultracompact, primordially-produced dark matter minihalos turn 
out to be especially promising in this respect. I illustrate how the results 
of these searches can be combined with constraints from accelerators and 
cosmology to produce a single global fit to all available data. Global fits 
in supersymmetry turn out to be quite technically demanding, even with 
the simplest predictive models and the addition of complementary data 
from a bevy of astronomical and terrestrial experiments; I show how ge- 
netic algorithms can help in overcoming these challenges. 

Key words: dark matter, supersymmetry, gamma rays, dwarf galaxies, 
stellar evolution, cosmological perturbations, phase transitions, statistical 
techniques 
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Preface 



This thesis deals with strategies for detecting and identifying dark mat- 
ter. Identification means somewhat more than simply establishing that 
dark matter is a particle with a certain mass, and belongs to one or an- 
other proposed phenomenological candidate class (WIMPs, axions, etc). 
Although this alone would be a tremendous achievement, it would be just 
the first step. In my opinion, it is the second that is actually the more 
interesting (and difficult) of the two. 

Rather, true identification means establishing which particle (or par- 
ticles) are responsible for dark matter, and characterising the actual field 
theory to which they belong. The papers included in this thesis all work 
towards achieving the first of these goals, but some (Paper IV and Paper 
V especially) also make contributions towards the second. 

Thesis plan 

This thesis is divided into two parts: Part I gives an introduction to the 
field and my work in particular, and Part II provides the included papers. 
Chapter 1 of Part I describes the observational evidence for dark matter, 
Chapter 2 details the various models proposed to explain it, and Chapter 3 
discusses the strategies which may be used to search for it. Chapter 4 gives 
a brief exposition of some of the most important uncertainties entering 
dark matter searches, and Chapter 5 summarises the results presented in 
the papers of Part II. 

The papers included in this thesis span quite a broad range of topics 
within the theme of dark matter searches. Many of them include introduc- 
tory sections with specific background material on the physics, existing 
literature and computational techniques relevant to their respective top- 
ics. Each paper also includes extensive discussion and interpretation of 
the results it presents; as one should expect, they are quite self-contained 
and largely self-explanatory. In most cases, it was I who was responsible 
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for actually writing those sections (or at least extensively revising them). 
Instead of simply repeating material already included in the papers them- 
selves, I have thus intentionally written the first four Chapters of Part I 
with the character of a textbook-level introductory review, so as to ex- 
plain and justify the broader scientific context in which the six papers 
have been written. In this undertaking I have not aimed for compre- 
hensive coverage of the literature. Full historical references and further 
details on the topics covered in Chapters 1-4 can be found in extensive 
review articles by Jungman, Kamionkowski & Griest (1996), Bergstrom 
(2000) and Bertone, Hooper & Silk (2005), as well as the recent reference 
volume edited by Bertone (2010). 

Unlike in some theses, the included papers of Part II should not be 
considered appendices. They are in fact the main body of this thesis, 
equivalent to Results, Discussion and Conclusion chapters in a traditional 
monograph. Their main results, figures and conclusions are summarised 
for the reader's convenience in Chapter 5, but the actual papers should 
be recognised as the definitive exposition of the work contained in this 
thesis. In short, the ethos I have adopted is that Part I should essentially 
contain information not provided already by Part II; there should be as 
little regurgitation as possible of content which is already contained in the 
six included manuscripts. Accordingly, Chapter 5 also places the results 
of Part II back into the broader context established by Chapters 1-4, and 
presents an outlook for future work. 

Contribution to papers 

My contribution to the included papers is as follows. For Paper I, I 
wrote and ran the stellar evolution code (as opposed to the static stellar 
structure code), and revised the manuscript. I discussed methods, results 
and conclusions together with the other authors. For Paper II, I designed 
most of the study, wrote the code and ran it using orbital and velocity data 
provided by Malcolm, and the results of neutrino simulations performed 
by Joakim. We discussed methods, results and conclusions together. I 
wrote the paper. I had the original idea for Paper III, designed the bulk of 
the study, wrote most of the paper and performed all calculations except 
for adiabatic contractions, which were done by Sofia. We interpreted 
results together. I designed almost all aspects of the work in Paper IV, 
with significant input on statistical matters from Jan. I wrote the software 
and ran it using data reduced by Christian, and wrote the paper. I 
interpreted results together with the other authors. For Paper V, I made 
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large contributions to the design and implementation of the study. I 
wrote sections of the paper and revised the remainder. I also organised 
the constituent physical codes into a consistent and functional starting 
point, from which Yashar added the genetic algorithm code and ran it. 
We discussed and interpreted results together with the other authors. 
For Paper VI, I helped design the study, wrote sections of the paper and 
computed TLUSTY model atmospheres. I discussed methods and results 
primarily with Erik, Fabio and Sofia. 

Pat Scott 

Stockholm, March 2010 
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Chapter 1 



The need for dark matter 



Dark matter has become such an estabhshed paradigm in modern astro- 
and particle physics that its existence is generahy accepted with httle 
explanation. Such blithe and widespread adoption of a cosmology in 
which an unknown matter plays the pivotal role has led some to become 
a little sceptical about its existence. It is worth reminding ourselves, 
from time to time, of the strong and compelling evidence on which this 
paradigm actually stands. As it so happens, we will also pick up a good 
deal of information along the way about what properties dark matter 
must have. 

1.1 Kinematics 

The earliest identification of dark matter came from the velocity disper- 
sions of galaxies within clusters. Zwicky (1933) noticed that the outer 
members of the Coma cluster were moving far too quickly to be merely 
tracing the gravitational potential of the visible cluster mass. The only 
way the observed velocities of the cluster members could be reconciled 
with the virial theorem was to postulate that the cluster also contained 
another large, but unseen, mass component: dark matter. Today, an 
equivalent technique is used to weigh clusters in large-scale X-ray surveys 
(e.g. Vikhlinin et al. 2006; Pacaud et al. 2007). In this case, because the 
hot. X-ray-emitting intracluster medium is entirely virialised, its temper- 
ature can be used to infer the total potential energy of the system, and 
therefore the cluster mass. Such surveys consistently confirm Zwicky's 
finding that clusters contain far more dark than luminous matter. 

It was not until similar effects were seen at the galactic level by Rubin 
and collaborators (Rubin & Ford 1970; Rubin et al. 1978, 1985) that the 



3 



4 



Chapter 1. The need for dark matter 




M33 rotation curve 



Figure 1.1. The observed rotation curve of the dwarf spiral galaxy M33, 
and the curve one would predict purely on the basis of the luminous disc. 
The fact that the curve is flatter than predicted suggests the presence of 
a halo of dark matter, extending to large galactic radii. From Bergstrom 
(2000). 



idea of dark matter really began to gain traction. Like the galaxies in 
the outer region of the Coma cluster, stars in the outer reaches of spiral 
galaxies are seen to rotate far more quickly than would be expected if each 
galaxy consisted only of the matter visible in stars and gas. An example 
of such anomalous circular velocities is shown in Fig. 1.1 for a dwarf spiral 
galaxy, M33. Here, the rotation curve expected purely on the basis of the 
stellar mass drops off markedly at long distances from the centre of the 
galaxy, whereas the observed curve is far flatter. This can be explained 
by the presence of an additional halo of dark matter, far more extended 
than the observed stellar disk. 

Similarly, the observed proportionality between the luminosity L of 
spiral galaxies and the maximum circular velocity V^max of their members. 



LcxVl,^, /3 = 3~4, (1.1) 

known as the Tully-Fisher relation (Tully & Fisher 1977), can be explained 
by dark matter and the virial theorem (Aaronson et al. 1979). The virial 
theorem leads to a direct correspondence between the total galactic mass 
M and T^aax- The self-similarity of hierarchically- formed dark matter 
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Figure 1.2. An image of a gravitationally lensed cluster (left), along 
with the corresponding mass map of the foreground cluster. In the image 
on the left, the foreground cluster can be seen in yellow, with the lensed 
image of the background cluster visible in blue. The reconstructed mass 
map shows a large, broad peak around the centre of the cluster which is 
not visible in the optical image, indicating the presence a massive dark 
halo. Credit: Greg Kochanski, Ian Dell'Antonio, and Tony Tyson (Bell 
Labs, Lucent Technologies). See also Tyson et al. (1998). 



halos (e.g. Navarro, Frenk & White 1996) and their resident baryonic 
disks causes this to be a power-law dependency. Assuming that for a 
given class of spiral galaxies (early-type, late-type, etc.) the mass-to- 
light ratio (M/L) is approximately constant, the relationship translates 
directly into a power-law dependency of L upon l^max- A similar argument 
holds for the analogous relation for elliptical galaxies, the Faber-Jackson 
relation (Faber &; Jackson 1976). 



1.2 Gravitational lensing 



According to General Relativity, the presence of any mass causes the 
space in its vicinity to curve. Curvature of space implies curvature of 
geodesies, which results in bending of light rays around massive bodies. 
In this way, light from background objects can be 'lensed' by massive fore- 
ground objects, as the light rays passing through the foreground objects' 
gravitational field are bent towards a distant observer. At the simplest 
level, this causes a background object to appear brighter than it otherwise 
would. This phenomenon provides a clean and unambiguous means by 
which to probe the mass distribution in a foreground object. 
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Figure 1.3. Composite image of the colliding 'Bullet Cluster', obtained 
with gravitational lensing. The lensing mass map is shown in blue, and 
the X-ray observations tracing the gas component are shown in pink. The 
offset of the mass and gas maps indicates that the mass of each cluster 
does not consist predominantly of gas, nor does it track the gas. Whilst 
the gas has collided with and exerted friction on the other cloud, leading 
to the characteristic 'punch-through' shape of the cluster to the right, 
the bulk of the mass has simply passed through the collision unimpeded. 
This proves that the majority of the matter in each cluster is collisionless 
dark matter. X-ray data: NASA/CXC/CfA/Markevitch (2006). Lensing 
Map: NASA/STScI; ESQ WFI; Magellan/U.Arizona/Clowe et al. (2006). 
Optical data: NASA/STScI; Magellan/U.Arizona/ Clowe et al. (2006). 



Distant clusters lensed by closer ones show evidence of extensive grav- 
itational lensing, far more than can be explained by the observed distribu- 
tion of luminous matter in the foreground cluster (e.g. Tyson et al. 1998; 
Massey et al. 2007). Fig. 1.2 gives an example of such a cluster, exhibiting 
a large, smooth mass component in the lensing map to the right, which 
is not accounted for by the distribution of luminous matter shown in the 
optical image to the left. 

Probably the most famous example of a dark matter lens is the so- 
called Bullet Cluster (Clowe et al. 2006), shown in Fig. 1.3. This object 
is in fact two clusters that have recently collided. In this case, not only 
does the lensing map (in blue) exhibit large amounts of dark matter not 
evident in the X-ray gas map (in pink), but the distribution of the two 
components also provides crucial insight into the properties of the dark 
matter. That is, the dark matter halos have passed straight through both 
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the gas clouds and each other, and appear essentially undisturbed after 
the collision. The gas clouds, on the other hand, have clearly exerted 
friction on each other as they have collided, leading to the characteristic 
ballistic shape of the rightmost cloud. This shows that dark matter does 
not necessarily have to track luminous matter in any way, and that it 
does not interact strongly with either gas or itself; this means that it is 
effectively collisionless. 

1.3 The large scale structure of the Universe 

On large scales, the Universe shows a wealth of structure: galaxies are 
gathered into clusters, clusters are part of superclusters, and superclusters 
are arranged into large-scale sheets, filaments and voids. This cosmic 
scaffolding has been revealed by large-scale surveys such as 2dFGRS (the 
2-degree Field Galaxy Redshift Survey; Colless et al. 2001) and SDSS (the 
Sloan Digital Sky Survey; Tegmark et al. 2004). Presumably, the pattern 
of galactic superstructure reflects the history of gravitational clustering of 
matter since the Big Bang. If dark matter were present during structure 
formation, it should have influenced the pattern of large-scale structure 
we see today. 

Large-scale cosmological 'N-body' simulations (e.g. Navarro et al. 1996; 
Springel et al. 2005; Diemand et al. 2007; Springel et al. 2008a; Diemand 
et al. 2008) demonstrate that the observed large-scale structure of lumi- 
nous matter could only have been formed in the presence of a substantial 
amount of dark matter. Furthermore, the bulk of dark matter must be 
both cold and non-dissipative for the correct structures to be produced. 
'Cold' in this context means that it moves very non-relativistically, and 
so has a short free-streaming length (less than the size of a gas cloud 
undergoing gravitational collapse, for example). Being cold means that 
the dark matter can gather gravitationally on small scales and so seed 
galaxy formation, but being non-dissipative prevents it from cooling and 
collapsing with the luminous matter, which would produce larger and 
more abundant galactic disks than are observed. Hot (highly relativistic) 
and warm (borderline relativistic) dark matter could still make up a frac- 
tion of the total dark matter, though just how large a fraction depends 
critically upon how warm that fraction is. 

Large-scale galaxy surveys not only provide information as to the 
amount and pattern of structure present in the Universe, but also a handle 
on the total mass contained within it (e.g. Tegmark et al. 2004; Cole et al. 
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2005; Percival et al. 2010). To a first approximation, the characteristic size 
of baryonic density perturbations which can survive until matter-radiation 
equahty (and therefore cohapse to form galaxies and clusters) is set by the 
mean density of baryons in the Universe (Silk 1968); smaller perturbations 
are destroyed by radiative and neutrino damping, larger ones fragment. 
In the presence of an additional matter component which gravitates but 
does not couple radiatively to the baryons, the power spectrum of per- 
turbations which ultimately survive to form galaxies is modified by the 
enhanced gravitational clustering. The resulting large-scale galaxy power 
spectrum is thus strongly dependent upon the total matter density of the 
Universe, and weakly dependent upon the fraction of matter contained in 
baryons. Recent surveys indicate a total matter (i.e. dark plus luminous 
matter) density of = Pm/Pc ^ 0.29 (Percival et al. 2010), where pc is 
the critical density required to close the Universe. 



1.4 Big bang nucleosynthesis 

A critical prediction of the hot Big Bang cosmology is that protons and 
neutrons were fused in the primordial fireball to create the light elements, 
as it cooled to temperatures of the order of an MeV. Modulo the effect of 
neutrinos, the resulting elemental abundances depend only on the nuclear 
reaction rates and the baryon-to-photon ratio (77) at the time. Labora- 
tory and theoretically-calculated nuclear reaction rates (e.g. Angulo et al. 
1999) can therefore be used to derive the primordial abundances of the 
elements as a function of 77 (e.g Coc et al. 2004; Steigman 2007; locco 
et al. 2009), as shown in Fig. 1.4. 

The parameter 77 is completely equivalent to rib^^ (up to a con- 
stant of proportionality), where fib is the baryon density of the Universe 
and h = i7o/100kms~^ Mpc~^ is the dimensionless Hubble parameter. 
Hq = v/d is the Hubble constant, describing the speed v at which galax- 
ies at a distance d appear to be receding due to the expansion of the 
Universe. With observations of the true primordial abundances of the 
elements and an independent measurement of /i. Big Bang Nucleosyn- 
thesis (BBN) therefore allows us to measure the primordial value of fib- 
Measurements of the Hubble parameter are abundant; the most widely- 
used comes from the Hubble Key Project {h = 0.72 zb 0.08; Freedman 
et al. 2001). The primordial abundances of the light isotopes are some- 
what more difficult to come by, requiring direct observations of extremely 
unevolved systems. 
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Because deuterium (D) is easily destroyed in stellar interiors and not 
produced in substantial amounts by other processes, its abundance shows 
a roughly steady decrease with time. The D/H ratio is often obtained by 
observing absorption lines in very old hydrogen clouds backlit by high- 
redshift quasars, and extrapolating backwards to nucleosynthesis. ^He, 
on the other hand, is created and destroyed in stellar interiors, the inter- 
stellar medium and the Earth's atmosphere, so its primordial abundance 
is very difficult to estimate reliably. The best upper limits to its primor- 
dial abundance come from emission line observations of the least evolved 
Galactic Hll regions, or the protosolar nebula. Similarly, ^He is produced 
ubiquitously in hydrogen fusion via the pp-chain in stars, so its abundance 
is substantially polluted relative to BBN; its higher abundance than ^He 
at least allows determinations using emission from extragalactic Hll re- 
gions. 

^Li is observed in absorption in the atmospheres of extremely metal- 
poor halo stars, where its abundance may eventually plateau with decreas- 
ing metallicity (Spite & Spite 1982). ^Li has also been seen in such stars 
(Asplund et al. 2006), though its presence can only be inferred from a very 
careful treatment of the isotopic broadening of spectral lines, along with 
the related three-dimensional effects of convection and departures from 
local thermodynamic equilibrium (LTE) on the line shapes. The heavier 
elements are rarely considered because their tiny primordial abundances 
make any reliable observational determination virtually impossible. 

The abundances of D and the He isotopes point very consistently to 
a baryon density of fib ^ 0.04. The corresponding prediction for the 
lithium isotopic abundances do exhibit some tension with observation: 
^Li is seen to be under-abundant relative to the BBN prediction, and 
^Li is not expected to have been produced in BBN at all. Both lithium 
isotopes are intrinsically unstable however, and subject to assorted pro- 
cessing in different stellar layers and the interstellar medium. The stellar 
determinations are also extremely difficult, depending crucially upon the 
input atomic data required for isotopic shifts and non-LTE calculations, 
and the fine points of detailed convective line shapes in stars known to be 
heavily effected by 3D corrections to their inferred abundances. Although 
it remains to be seen if the discrepancies will indeed be borne out by fu- 
ture observations, there is speculation that they might be explained by 
non-standard BBN caused by particle physics beyond the standard model 
(BSM; beyond the SM), including some specific dark matter models (e.g. 
Jedamzik & Pospelov 2009). 
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Figure 1.4. Abundances of the light isotopes produced in Big Bang 
nucleosynthesis (BBN), as a function of the baryon-to-photon ratio 7710 = 
10^°?7, which directly corresponds to the baryon fraction of the universe 
comparing such predictions with observations of very 
old astronomical systems, one can use BBN to derive Qb- Abundances 
are given as mass fractions. The bands between the two curves for each 
isotope indicate the degree of uncertainty arising from nuclear reaction 
rates and other theoretical sources. Yp refers to primordial ^He. From 
Steigman (2007). 



The independent measurements of fib from BBN and from large- 
scale structure together provide incontrovertible evidence for the exis- 
tence of dark matter. Since all luminous matter essentially consists of 
baryons, fim ^ 0.29 and ^ 0.04 together imply that the remaining 
fi leftover ^0.25 must be dark matter. Furthermore, we have another in- 
valuable piece of information about its nature: we see immediately that 
dark matter must be predominantly non-baryonic. Being non-baryonic 
allows the possibility that the dark matter is not capable of interacting 
with photons electromagnetically (something baryons clearly are capable 
of), which sits very well with the fact that it does indeed appear dark. 
This is also consistent with dark matter being dissipationless; were it able 
to absorb and re-emit photons, it would obscure stars within distant ha- 
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los, as well as radiate away angular momentum and collapse with baryons 
to form stellar and galactic disks. 

It is worth mentioning that some amount of baryonic dark matter also 
remains unaccounted for today in the figure fib ^ 0.04, in that it has not 
been directly observed in surveys of gas or galaxies. Nevertheless, we 
know this to be far less than what exists in non-baryonic dark matter. 
For the purposes of this thesis, 'dark matter' is taken as shorthand for 
the dominant, non-baryonic component. 

1.5 Concordance cosmology 

In a sense, the evidence presented in Sects. 1.1-1.4 is already sufficient to 
conclude that dark matter definitely exists, and must be neither baryonic 
nor able to interact substantially with baryons, photons or itself. Final 
confirmation of this comes from observations of the cosmic microwave 
background (CMB). 

The large-scale distribution of galaxies shows density variations on 
scales which reflect the situation at matter-radiation equality. Similarly, 
the temperature inhomogeneities of the CMB exhibit characteristic scales 
which reflect the situation shortly after, at the time of recombination (see 
e.g. Hu & Dodelson 2002; Samtleben et al. 2007). Unlike the galaxy power 
spectrum however, the angular temperature anisotropics of the CMB ex- 
hibit the primordial density perturbations of the coupled baryon-photon 
fluid directly; following recombination, the angular distribution of the 
massless photons is held nearly uniform by radiation pressure, whereas 
the matter collapses due to gravity. The angular power spectrum of the 
CMB observed today is therefore sensitive to the full range of cosmolog- 
ical parameters which play a role in the evolution of the baryon-photon 
oscillations: the total energy density of the Universe, the baryon fraction 
and the spectral shape of the primordial perturbations. It is also sensi- 
tive to the large-scale geometry of the Universe since recombination, as 
the observed angular diameter of the characteristic scales frozen into the 
CMB at recombination depends upon the geometry of the space through 
which they have travelled to reach us. 

Thanks to the extremely high resolution of recent CMB missions, 
fits to the microwave background provide accurate measurements of the 
matter density of the Universe, as well as the baryon fraction. The 
7-year Wilkinson Microwave Background Probe (WMAP) results (Lar- 
son et al. 2011) give posterior mean values of = 0.267 ± 0.026, 
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Figure 1.5. The current concordance model of cosmology, indicating 
that the Universe is dominated by dark energy and dark matter, and es- 
sentially flat. These figures illustrate the agreement between cosmological 
fits to BBN, observations of the cosmic microwave background (CMB), 
large scale structure (BA0/2dFGRS) and Type la supernovae. The fit 
is so good that any one of the data sets could be removed and essentially 
the same cosmology would be inferred from the remainder. Contours in 
the figure on the left (a) indicate 1-, 2- and 3cr plausible ranges, whilst 
areas plotted on the right (h) are la ranges. Left figure from Kowalski 
et al. (2008, reproduced by permission of the AAS), right figure adapted 
from Tegmark et al. (2001) with additional 2dF results from Cole et al. 
(2005). 
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Jib = 0.0449 ± 0.0028, I^dm = 0.222 zb 0.026. These results are in excel- 
lent agreement with those obtained from BBN and large-scale structure, 
confirming the need for non-baryonic dark matter beyond any doubt. 

The CMB fits also indicate that the Universe is approximately flat, 
requiring some form of 'dark energy' to make up the remaining Vt^ = 
0.734 zb 0.029. Although difficult to explain theoretically, the need for 
this component was confirmed by observations that distant supernovae 
are apparently dimming, indicating that the expansion of the Universe is 
accelerating (Riess et al. 1998; Perlmutter et al. 1999). The CMB, Type 
la supernovae, large scale structure and BBN together paint an entirely 
self-consistent picture of the Universe in which we live: baryons make 
up little more than a few percent of the energy budget of the Universe, 
non-baryonic dark matter about a quarter, and dark energy the remain- 
ing three quarters. The remarkable level of consistency between these 
different data sets can be seen in Fig. 1.5; were any one of the data sets 
to be removed, the model would stand perfectly well constrained with the 
remaining three. A combined fit to all data (Komatsu et al. 2009) gives 

= 0.274 ± 0.013, J^b = 0.0456 ± 0.0015, J^dm = 0.228 ± 0.013 and 

= 0.726 ±0.015. 
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As we have seen in the previous chapter, the identity of dark matter 
is not entirely unconstrained. Firstly, it must be massive, because we 
see its gravitational influence. We know that it must of course also be 
dark, which implies that it can neither absorb nor emit photons at rates 
comparable to normal matter. This means that it must have a very 
small effective coupling to photons, either because it exists in some very 
special geometry that prevents it coming into contact with them very 
often, or due to its intrinsic particle properties. In this case, it must be 
electrically neutral or possess some extremely small fractional charge, as 
any particle with charge q couples to photons via tree-level QED 
processes, leading to cross-sections proportional to q'^a^ (where a is the 
fine-structure constant). 

We know that dark matter must be non-dissipative in order to pro- 
duce the correct structures in the early Universe. This also implies that 
there can be no substantial (i.e. electromagnetic-strength or stronger) 
interaction between dark matter and normal matter. If there were, en- 
ergy would be efficiently transferred from the dark matter to the normal 
matter and radiated away as photons, allowing the dark matter to form 
disks along with normal matter. Dark matter is thus effectively collision- 
less with respect to all normal matter: baryons, leptons and photons. We 
know then that it must be non-baryonic, because it does not couple to 
photons or baryons; this is spectacularly confirmed by the anisotropics 
in the CMB. Observations of systems such as the Bullet Cluster provide 
even more direct proof that dark matter does not couple substantially to 
baryons, and even show that dark matter must be effectively coUisionless 
with respect to itself. 
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Dark matter must also be cold in order for structure formation to 
proceed correctly, which places limits on the mass(es) of its constituent 
particle(s) in some cases. Finally, any dark matter candidate must be 
produced with the correct relic abundance to plausibly explain observa- 
tions that f^DM = 0.23 today. This also means that it must be stable, or 
at least very long-lived, in order to have persisted in significant numbers 
to the present day. 

2.1 The good, the bad and the ugly 

A great many theories have been put forward for the identity of dark 
matter. A good candidate theory fulfils all the requirements above, fits 
all other experimental constraints (most of which are described below or 
in Chapter 3), requires minimal arbitrary choices in its parameter values, 
and makes some testable predictions. The last two points are debatable; 
one is the requirement that the candidate presents a subjectively 'natural' 
solution to the problem, and the other simply makes the theory practically 
useful. Mainly due to the non-baryonic requirement, all theories which 
qualify as 'good' today fall into the realm of particle dark matter: namely, 
that dark matter consists of some as-yet undiscovered particle. The most 
notable theories of this kind are outlined in Sects. 2.2-2.7. 

A theory that was rather popular in astronomical circles for a time 
was that of MACHOs (Paczynski 1986), MAssive Compact Halo Objects. 
These were hypothesised to be dim but otherwise quite normal astrophysi- 
cal objects, which populated the outer halos of galaxies. Known candidate 
objects were red, white and brown dwarfs, black holes and substellar-mass 
objects (in effect, hostless planets). These systems of condensed baryons 
would appear dark because they would be so dim, but would also be ef- 
fectively transparent to background light because their compact nature 
would cause the effective cross-section for light absorption per unit mass 
across an entire halo to be very small. The obvious trouble with MA- 
CHOs was that they would have been baryonic dark matter, violating the 
non-baryonic requirement. There was also never any convincing argument 
for why such objects should preferentially congregate in the outer halo, 
whereas other more luminous stars should not. Nonetheless, MACHOs 
were popular because they were testable (not because they were logical): 
microlensing searches towards the Magellanic Clouds (Alcock et al. 2000; 
Tisserand et al. 2007; Wyrzykowski et al. 2009) promptly ruled MACHOs 
out as the dominant contributor to dark matter some years ago. 
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Primordial black holes (PBHs) are another dark matter candidate 
with significant problems. These would be formed from small-scale pri- 
mordial density perturbations which are so strong that they cause the 
entire horizon mass at the time (or at least a very substantial fraction 
of it) to collapse directly into a black hole. The fluctuations would be 
produced before matter-radiation equality, either during inflation or by 
phase transitions in the early Universe. Because PBHs form before BBN 
or the CMB, the baryons of which they consist are effectively sequestered 
from the rest of the Universe, allowing them to act like non-baryonic dark 
matter. Like MACHOs, they would escape detection via shadowing ef- 
fects today simply by virtue of their compactness. The difficulty with this 
scenario is that the density contrast 

g I /^perturbation Pambient | ^2 

Pambient 

required to form a PBH is of the order of 5 > 30% (e.g. Green et al. 
2004b). In comparison, the initial density perturbations from inflation 
were about S ^ 10~^. A very bottom-heavy spectrum of perturbations is 
then required in order to provide enough power on small scales to produce 
a substantial number of PBHs without violating the level of large-scale 
anisotropy seen in the CMB. 

Another problem is that PBHs can be produced with virtually any 
primordial abundance, creating a major fine-tuning problem. Their relic 
density is obtained by integrating the probability distribution of primor- 
dial density perturbations between their formation threshold {6 ^ 30%) 
and 5=1. Assuming e.g. a Gaussian power spectrum of perturbations for 
simplicity (e.g. Green & Liddle 1997), this produces a density at matter- 
radiation equality of 

f^PBH(MH) = r , exp f- ,.,f, d5. (2.2) 

^0.3 



2^a(MH) V 2a(MH 

Here (j(Mh)^ is the variance of perturbations at the time when the horizon 
mass is Mh- With a scale-independent perturbation spectrum of index n. 



^(Mh) oc M^""'^. (2.3) 

On CMB scales, n is close to 1 (Komatsu et al. 2011), but at small scales, 
it could be very different. The resultant relic density from Eq. 2.2 is ex- 
tremely sensitive to n, even for a single epoch of formation (i.e. a single 
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value of Mh). The same is essentially true for other reasonable choices 
of the power spectrum and scale-dependence. See Paper III for a more 
detailed discussion. To the best of my knowledge, there is no known 
mechanism which naturally produces exactly the values of n and Mh 
required to reproduce the observed cosmological density of dark matter 
with PBHs. Even more disastrously though, PBHs are mostly disfavoured 
now by experiment. More massive PBHs should have been picked up by 
microlensing searches, and the Hawking radiation produced in the evapo- 
ration of low mass PBHs should have created effects on the extragalactic 
diffuse gamma-ray background, BBN, the CMB or cosmic ray data that 
have not been seen (see Josan et al. 2009, for an extensive compilation of 
limits). 

SM neutrinos were also considered viable dark matter candidates for 
a time (e.g. Steigman et al. 1978), since they are massive, stable, non- 
baryonic and do not interact electromagnetically with other matter. With 
the current limits on their masses however (e.g. Schwetz et al. 2008), they 
are far too light to meet the definition of cold dark matter, and can only 
have a quite small cosmological abundance {0.^ < 0.03; Bergstrom 2009). 

For a time, the notion of a heavy fourth-generation SM-like neutrino 
was popular as a candidate for dark matter. In order to achieve the correct 
relic abundance, its mass would have to be either 1-50 eV, or more than 
about a few GeV (Jungman et al. 1996; Bergstrom 2000); the former is 
ruled out for SM-like neutrinos by the Tremaine-Gunn limit (Tremaine & 
Gunn 1979), and the latter is excluded by direct searches for dark matter 
(Sect. 3.1). 

To some, the MOND paradigm (MOdified Newtonian Dynamics; Mil- 
grom 1983; Bekenstein 2004) is appealing. This approach involves modi- 
fying gravity in an ad hoc fashion at long distances in order to fit galactic 
rotation curves without the need for dark matter. In this sense, one could 
say that it is quite a successful empirical fitting formula for parameteris- 
ing the effects of dark matter on the kinematic observables discussed in 
Sect. 1.1. As a physical theory however, it has no basis to speak of, and 
fails entirely to explain any of the additional evidence for dark matter 
presented in Sects. 1.2-1.5. 

2.2 Axions 

Because gluons should be pure gauge fields at spatial infinity, the QCD 
vacuum possesses a rather complex structure (Peccei 2008). Being a pure 
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gauge field implies only that these boundary values must be drawn from 
the set of field configurations which can be obtained by gauge transfor- 
mations on the null field (i.e. 0). The vacuum structure arises because 
the pure gauge boundary condition introduces a freedom in the choice 
of boundary field, which translates into similar freedom in the choice of 
QCD vacuum. Because of this vacuum structure, the QCD Lagrangian 
picks up an additional effective term 



tensor for the ath gluon, and Ga/^u = 2^^uoL^Ga is its dual. Unlike the 
rest of the QCD Lagrangian, Cq does not conserve CP. This is a problem, 
as the strong interaction is known empirically to conserve CP rather well. 
This additional term would induce effects which have not been observed, 
like a non-zero electric dipole moment for the neutron. The limits on 
such a moment to date constrain the vacuum angle 9 to be less than 10~^ 
(Peccei 2008). Understanding why 9 is so small is known as the 'strong 
CP problem'. 

One widely held suspicion is that the solution to the strong CP prob- 
lem lies in the existence of an additional spontaneously-broken, global 
chiral U{1) symmetry of the SM Lagrangian, known as the Peccei-Quinn 
symmetry (Peccei & Quinn 1977). The axion is the Goldstone boson of 
this broken symmetry (see e.g. Kim & Carosi 2010, for a review). It pos- 
sesses a potential with a minimum that naturally sets the field to a value 
that cancels Cq. This is true for essentially any value of the Peccei-Quinn 
symmetry breaking scale relieving the naturalness problem posed by 
the smallness of 9. 

The low-energy phenomenology of axions is set almost entirely by their 
decay constant fa — Va/N. Here N is an integer which gives the degree of 
the colour chiral anomaly of the Peccei-Quinn symmetry (a corresponding 
electromagnetic chiral anomaly also exists, and is typically indexed with 
another integer E). The axion mass for example is given by (Sikivie 2008) 





(2.5) 



Axions possess a vertex with two photons (Raffelt 2008), 




(2.6) 
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where z = m^/m^i is the ratio of up and down quark masses, E and 
B are the electric and magnetic fields, and a is the axion field itself. 
The two photon vertex not only allows axion decay to two photons, but 
allows axion conversion to photons (and vice versa) in the presence of 
electromagnetic fields. This feature is used as the prime phenomenological 
means for searching for axions, as discussed in Sect. 3.6. 

Certain classes of axions, dubbed 'invisible axions', make very good 
dark matter candidates because they interact extremely weakly with nor- 
mal matter. This class of axion involves a very high Peccei-Quinn break- 
ing scale, and therefore very light axions. Such axions have virtually 
no kinematically-accessible decay channels, so are stable on cosmologi- 
cal timescales. These axions constitute cold dark matter despite being so 
light, because they are never in thermal equilibrium in the early Universe, 
so are never heated to relativistic temperatures along with the other forms 
of matter. 

There are two commonly discussed mechanisms for the primordial 
production of dark matter axions: vacuum misalignment and string decay. 
Depending on the timing of Peccei-Quinn breaking relative to inflation, 
the correct relic density of dark matter axions can be obtained either only 
for a very specific value of fa (^7 x IQ-*^^ GeV), or for a tightly correlated 
corridor of values in the fa-^\ plane (Visinelli & Gondolo 2009). The 
variable Lp\ is the initial value of the misalignment angle between the 
axion field and its minimum; in this model, it is the initial offset from 
the minimum which causes the field to oscillate and create particles. In 
this sense axions have some fine-tuning issues of their own as dark matter 
candidates, though they are not nearly as severe as for e.g. primordial 
black holes. 

2.3 WIMPs 

The most widely-studied dark matter candidates are the Weakly Inter- 
acting Massive Particles (WIMPs). WIMPs interact with SM fields only 
via the weak nuclear force, making them non-baryonic and electrically 
neutral by definition. WIMPs must carry some sort of conserved quan- 
tum number to keep them stable on cosmological timescales. This is 
usually achieved by making the WIMP the lightest member of a mat- 
ter sector which is charged under some discrete symmetry. More often 
than not, this symmetry needs to be essentially imposed by hand on the 
underlying theory. Being members of SU{2)i^ multiplets, in the absence 



2.3. WIMPs 



21 



of any other symmetries which might force their masses down, WIMPs 
should naturahy acquire masses within a few orders of magnitude 
of the SU{2)i^ X U{1)y electroweak symmetry-breaking scale (i.e. a few 
GeV or TeV). This makes them sufficiently heavy to constitute cold dark 
matter even if they have been produced thermally in the early Universe. 
Examples of WIMPs include the lightest neutralino in supersymmetry 
(Sect. 2.4.2), the lightest Kaluza-Klein (KK) particle (Servant & Tait 
2003) and an additional inert Higgs boson (Barbieri et al. 2006). 

2.3.1 Thermal relics 

A thermal relic is a particle whose cosmological abundance is set by ther- 
mal production of the particle in the early Universe. 'Thermal production' 
actually refers to chemical decoupling of a species from the primordial 
particle soup created in the Big Bang. Here 'chemical' has the mean- 
ing of particle creation and destruction by quantum processes, not the 
usual meaning of creation and destruction of molecules and free atoms by 
atomic interactions. True thermal decoupling on the other hand, other- 
wise known as kinetic decoupling, refers to the time when the velocities 
of a relic species cease to reflect the temperature of the Universe. 

Chemical decoupling, or 'freeze out', happens when the expansion of 
space eventually overcomes the rate at which the species is created and 
annihilated in interactions with other particles. The primordial fireball 
necessarily begins in chemical (and thermal) equilibrium at some very 
high temperature T^. The populations of different particle species are set 
by the equilibrium rates of particle creation and annihilation at Ti, and 
their velocities follow a Maxwellian distribution with temperature Ti. As 
the Universe expands, it cools, and for a time, the equilibrium populations 
of the particles adjust accordingly as chemical equilibrium is maintained. 
At some point in the expansion, when T — T^^ the creation and annihi- 
lation of a particular species will cease to be able to keep pace with the 
expansion, and the particle will fall out of chemical equilibrium. If it is 
stable, at this point the comoving density of the species will become essen- 
tially fixed. ^ Similarly, when elastic collisions between particles become 
sufficiently infrequent, both due to the increasing scale and decreasing 
temperature of the Universe, it will fall out of kinetic equilibrium. 

^This is analogous to the transition from LTE to non-LTE populations of atomic 
energy levels in astrophysical plasmas as one considers progressively lower gas densities 
or temperatures. 
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The time of chemical freeze out therefore sets the final population of 
any stable relic, and the time of kinetic decoupling sets its temperature. In 
general, kinetic decoupling of WIMPs happens after chemical decoupling 
(see e.g. Bringmann 2009). With weak-scale masses and creation/anni- 
hilation cross-sections, ^ 0.05m^ (Jungman et al. 1996). As kinetic 
decoupling happens at comparable or even lower temperatures than this, 
WIMPs freeze out at sub-relativistic energies. They are therefore guaran- 
teed to be moving non-relativistically at the time of structure formation, 
so qualify as cold dark matter. 

The resultant abundance of the relic species depends upon the his- 
tory of the expansion rate up until freeze out, and the rate at which the 
species annihilates into lighter particles. The expansion rate is H = a/a 
(where a is the scale factor of the Universe and the dot denotes time 
differentiation), and is obtained from the cosmological model. To a first 
approximation, the evolution of the population is thus described by the 
Boltzmann equation (Jungman et al. 1996) 

^ + 3Hn^ = {av) (n^, ^q' - n^^) (2.7) 

where n^^eq and are the equilibrium and actual number densities of 
the species and (av) is the thermally-averaged product of the relative 
velocity and the total annihilation cross-section for XX ~^ other particles. 
Eq. 2.7 holds whether x is Majorana or Dirac, though in the Dirac case it 
is only valid if there is no particle-antiparticle asymmetry to begin with, 
meaning that the total particle population is given by 2n^. 

Eq. 2.7 can be solved numerically in order to obtain the relic density 
of any species of interest; Fig. 2.1 gives solutions for the chemical freeze 
out of some example WIMPs. To a good approximation (Jungman et al. 
1996; Bergstrom 2000), for particles with masses in the relevant range for 
WIMPs 

^ (av) ^ ^ 

As seen in Eq. 2.8 and Fig. 2.1, larger self-annihilation cross-sections cause 
the particle to freeze out later, and so exhibit a lower relic density. In 
practice, co-annihilations between the WIMP and other particles nearly 
degenerate with it in mass, as well as effects such as annihilation thresh- 
olds and resonances, can make the calculation a little bit trickier than 
Eqs. 2.7 and 2.8 might suggest. Dedicated computer packages have been 
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Figure 2.1. Freeze out of a stable thermal relic. At early times, the 
relic particle is in chemical and thermal equilibrium with other parti- 
cles, with its abundance set by the equilibrium rates of particle creation 
and annihilation. Because of thermal coupling, these are set entirely by 
the temperature of the cooling Universe. At later times, the expanding 
Universe cools to the point where the equilibrium density is too low to 
maintain chemical contact with other particles. At this point, the particle 
falls out of chemical equilibrium with the rest, and its comoving density 
becomes fixed. The density at which this occurs depends upon the annihi- 
lation cross-section of the particle; larger cross-sections allow equilibrium 
to be maintained for longer, resulting in lower relic abundances. From 
Jungman et al. (1996), modelled on a figure from Kolb & Turner (1990). 



developed to solve Eq. 2.7 in detail, for a variety of specific WIMP models 
(e.g. Gondolo et al. 2004; Belanger et al. 2007). 

Eq. 2.8 shows that the canonical annihilation cross-section implied 
by the cosmological observations discussed in Chapter 1 is {av) ^ 3 x 
10~^^cm^s~^. The natural value of the annihilation cross-section for 
a weakly-interacting particle can be estimated as {av) ^ ^^/^weak^ ^ 
10"^GeV"2 = 10"25cm^s-\ where £;weak ^ 100 GeV is the electroweak 
scale. ^ Here we have an amazing coincidence: purely by postulating that 



^That this is approximately true can be verified by e.g. considering Z- mediated 
X self-annihilation, such that \M\^ ^ g^E^^ /niz^ , and using Eq. 4.84 of Peskin & 
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dark matter is a stable weakly-interacting particle, we have predicted the 
relic density to within an order of magnitude. The naturalness with which 
WIMPs can provide the correct relic density has lead this to be dubbed 
'the WIMP miracle'. Whilst far from definitive, it is a very strong hint 
that WIMPs might be the correct solution to the dark matter puzzle. 



2.3.2 The hierarchy problem 

Indeed, the process of chemical freeze out is so generic that if a stable 
neutral particle exists around the electroweak scale, then we really should 
expect it to be dark matter. But what reason do we actually have for 
believing that such stable particles might exist around or just above the 
electroweak scale? This is a question of why we think there might be 
new physics at such an energy scale, as new physics generically means 
new particles (though not necessarily stable ones). There are a number 
of reasons to expect that new physics exists beyond the SM, since the 
SM cannot explain phenomena such as the net baryon asymmetry of the 
Universe, neutrino masses, the imperfect unification of the known gauge 
couplings, and the spectrum of particle masses - but the one compelling 
reason to believe that new physics exists specifically at the TeV scale is 
the mass of the Higgs boson. 

As a fundamental scalar, the mass of the Higgs is not protected by any 
symmetry; quadratic divergences generated by fermionic loop diagrams 
such as Fig. 2.2 have no apparent cancellation or suppression, nor can 
they be absorbed into the renormalisation group running of any coupling 
constant. For fermions, the case is different: similar mass corrections are 
controlled by the approximate chiral symmetry of the fermionic part of 
the SM Lagrangian (Luty 2005). For the gauge bosons, gauge symmetry 
suppresses the divergences and shifts them into the running of the gauge 
couplings. For the SM Higgs however, the quadratic divergences remain 
bare, causing extreme corrections to its mass. There is therefore no reason 
to believe that the SM Higgs should have a mass any lower than the GUT 
scale. That is, except for the pesky detail of its fundamental role in 



Schroeder (1995) to arrive at 

(2.9) 



where g is the electroweak couphng and G is the Fermi coupUng constant. In the case 
where the incoming xs are moving non-relativistically (as WIMPs do), E^ ^ ^ 
E^ea,]^^ giving the zero- velocity limit {(tv)q ^ 10~^ GeV~^ = 10~^^ cm^ s~^. 
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Figure 2.2. An example one-loop contribution to the Higgs mass from 
a fermionic loop. The amplitude of this diagram is quadrat ically diver- 
gent, so the Higgs mass is sensitive to some large mass scale (which may 
be expressed as the cut-off used to regularise the diagram's amplitude, 
amongst other presentations; see e.g. Luty 2005). 



electroweak symmetry breaking and the generation of mass for the entire 
spectrum of SM fermions, which it needs an elect roweak-scale vacuum 
expectation value (VEV) to achieve. How then is the Higgs mass stabilised 
at the weak scale? This problem is known as the gauge hierarchy problem, 
and strongly suggests that some new physics must kick in at energies just 
above the electroweak scale, allowing us to calculate the effective Higgs 
mass below this scale by regularising the amplitudes of the divergent 
diagrams with a cut-off at A ^ 1 TeV. 

The hierarchy problem is a very generic motivation for new physics 
at or just above the electroweak scale. That is to say, it is independent 
of what the nature of this new physics might be, so there is no a priori 
reason for us to naively expect that it necessarily has anything to do with 
dark matter. The fact that its resolution would appear to require new 
physics at exactly the right scale for the associated particles to be WIMPs 
is thus another amazing coincidence, and a second very strong hint that 
WIMPs constitute the majority of dark matter. 

2.4 Supersymmetry 

Although this chapter is primarily intended as something of a bestiary 
of dark matter models, here it is worth detouring slightly into the the- 
ory of supersymmetry (SUSY). This is because SUSY not only provides 
multiple examples of WIMPs, but because it is also exactly the sort of 
'new physics' discussed in the previous section as necessary to solve the 
hierarchy problem. 
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2.4.1 General features 

The basic idea of SUSY is to give every SM particle one (or more) super- 
partner (s), which differ from their SM counterparts only by their spins 
(see e.g. Martin 1997; Baer & Tata 2006; Aitchison 2007, for reviews; 
the bulk of the treatment in this section follows Baer & Tata). Every 
SM fermion has a sfermionic superpartner (which is a boson) and every 
SM boson has a bosino superpartner (which is a fermion). The symmetry 
is a continuous extension of the Poincare group to fermionic (i.e. anti- 
commuting) degrees of freedom. Supersymmetry is thus essentially an 
extension of special relativity into a four-dimensional fermionic 'super- 
space'; the fermionic dimensions can be viewed either as four additional 
ant i- commuting dimensions, or simply as fermionic extensions of the ex- 
isting four. Correspondingly, the generators of SUSY transformations are 
themselves spinorial, in contrast to the vector and scalar generators asso- 
ciated with traditional spacetime and internal symmetries. It is entirely 
possible for there to exist multiple distinct supersymmetry generators, 
and therefore multiple supersymmetries (known as > 1 supersymme- 
tries, where N is the number of generators). = 1 supersymmetry is 
the only version which permits chiral fermions, so known phenomenology 
dictates that supersymmetry can only be of the A^ = 1 variety at low 
energies. 

It is often said that SUSY provides a non-trivial way of unifying space- 
time and internal symmetries, since the anti-commutation relation for the 
spinor components of the SUSY generator Q 



gives the generator of spacetime translations, P^.^ This means that the 
combination of any two SUSY transformations is equivalent to a spacetime 
translation. So, the supersymmetry algebra mixes spacetime transforma- 
tions with boson- fermion exchange, which would seem to be something 
of an 'internal' operation because it directly modifies particle proper- 
ties. This interpretation seems unreasonable though, as the swapping of 
bosons for fermions or vice versa involves modifying only the particle's 
spin, which is arguably no less of a spacetime property than momentum 
or position. Even if supersymmetry were gauged (more on this later), the 

^Here I am working with 4-component spinor notation as per Baer & Tata (2006). 
The bar denotes the spinor adjoint Q = Q^7°, and 7^ are the standard Dirac gamma 
matrices, forming a representation of the Chfford algebra of Minkowski space. 



{Qa,Qb} = {lt.)abP'' 



(2.10) 
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Figure 2.3. An example one-loop contribution to the Higgs mass from a 
sfermionic loop. The amplitude of this diagram has exactly the same sen- 
sitivity to the high-scale cut-off as the fermionic loop diagram of Fig. 2.2, 
except for a factor of — |. In unbroken supersymmetry, two of these di- 
agrams appear for each fermionic one, exactly cancelling the fermionic 
loop's contribution to the Higgs mass. The same occurs to all orders 
in perturbation theory. This removes the sensitivity to the high scale, 
stabilising the Higgs mass against higher-order corrections. 



resulting gauge theory of gravity would still seem to be an entirely space- 
time symmetry, and all the other gauge symmetries of particle physics 
still entirely internal ones. This is effectively what was said by the SUSY 
version (Haag et al. 1975) of the celebrated Coleman-Mandula 'no-go' 
theorem (Coleman & Mandula 1967): it is not possible to have a larger 
spacetime symmetry than the super-Poincare algebra, so there are no 
non-trivial ways to unify spacetime and internal gauge symmetries into 
the same group. 

SUSY solves the hierarchy problem in a very simple and elegant way: 
because every fermion has a bosonic superpartner and vice versa, every 
divergent loop diagram containing a SM fermion like Fig. 2.2 is matched 
by corresponding scalar sfermion loop diagrams like Fig. 2.3. In this 
case there are actually two scalar loop diagrams, because chiral fermions 
have two degrees of freedom, so their scalar superpartner fields must be 
complex, resulting in two oppositely charged sfermions for each chiral 
fermion. Because closed fermionic loops introduce a factor of —1 to their 
diagram's amplitude (whereas bosonic loops do not), and the particles 
are otherwise identical, the (regularised) contributions of the two copies 
of the scalar loop diagram cancel that of the fermionic one. The same 
is so to all orders in perturbation theory, so the Higgs mass receives no 
extreme corrections, and is stabilised somewhere near the SUSY-breaking 
scale. 

SUSY must clearly be a broken symmetry, as we have not observed any 
of the superpartners to date (known particles cannot be the superpartners 
of each other, as they do not fit correctly into the supermultiplet structure 
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required by the supersymmetry algebra). Except for some very specific 
scenarios where particular superpartners happen to be especially sneaky, 
we can infer that they must all be considerably more massive than the 
SM particles, as they have yet to show up in e.g. collider experiments. We 
would expect then that the SUSY particles should generally have masses 
of around the SUSY-breaking scale. The cancellation of divergences thus 
cannot be exact following SUSY-breaking, as particles and sparticles must 
have different masses. We can then conclude that if SUSY is a viable 
solution to the hierarchy problem, it must be broken at about a TeV, 
otherwise the cancellation of the different loop contributions would again 
require an unexplained fine-tuning. In fact, another appealing feature 
of SUSY is that it can actually cause electroweak symmetry-breaking 
(EWSB), in that the Higgs mass term in the SU{2)i^ x U{1)y Lagrangian 
can be driven negative by SUSY renormalisation group evolution - but 
only if the SUSY-breaking scale is about a TeV. So, we might expect 
EWSB to have been induced by SUSY breaking, and hence taken place 
at only a slightly lower energy than SUSY-breaking itself. Since we know 
the electroweak scale to be a few hundred GeV, this gives further credence 
to the idea that SUSY-breaking should occur at around a TeV. 

Actually, the cancellation of quadratic divergences depends a lot more 
crucially upon the gauge coupling constants in the SM and SUSY sector 
being identical than the particle masses. It is possible to introduce explicit 
mass terms to the SUSY Lagrangian which break SUSY 'softly', in that 
they give the sparticles different masses to the SM particles, but do not 
introduce any quadratic divergences. 

The simplest prescription then for writing down a minimal supersym- 
metric theory is to take the supersymmetrised version of the SM La- 
grangian, >CsuSY-SM5 and augment it with all possible renormalisable soft 
SUSY-breaking terms (Baer & Tata 2006) 




(2.11a) 



(2.11b) 



(2.11d) 



(2.11c) 
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+ 



+ 



+ {Ae)ijL'}Hdael. + h.c^ (2.11e) 

+ (Ce)i,L^i/l4. + h.c.]}. (2.11f) 

Here the generation indices i and j are explicitly summed over, whilst 
summation is implied over the SU{2)i^ indices a^h — 1,2 and the gauge 
generator indices A = 1..3 and B = 1..8. Tilded operators denote 
superparticle fields (Qj, u\^j^ etc.), where S^, Wa and qb are the super- 
partners of the familiar SM gauge bosons. Other capitalised fields are 
SU(2)i^ doublets, whilst lowercase fields are singlets: Qi = (uu^du)^ ^ 
Li = {eu,i^e.L)^, Hu = {H+,Hy^ and = {HJ,H^)^. The dou- 
blet definitions for the superpartner fields are simply tilded versions of 
the same. Notice that there are now two Higgs doublets instead of the 
normal one from the SM. After EWSB, this gives rise to two additional 
charged Higgs bosons and a further two neutral Higgs compared the SM 
(recall that three real degrees of freedom are lost to the gauge bosons 
in EWSB, so two doublets means 8 — 3 = 5 real degrees remaining, in- 
stead of 4 — 3 = 1 real degree in the SM). The two doublets give masses 
separately to the up-type and down-type quarks, hence the subscripts u 
and d] the up-type doublet corresponds to the SM Higgs. The down-type 
doublet is required because the up-type Higgs cannot couple to the down- 
type quarks without breaking the SUSY-invariance of the action. Unless 
multiple Higgs doublets are present, Higgsinos also ruin certain anomaly 
cancellations which exist in the SM. 

The first two terms in Eq. 2.11 (2.11a and 2.11b) give explicit masses 
to the gauginos via the real parameters Mi, M2, M3 and M{, M2, M3. The 
second of these terms violates CP, so M[,M2 and M3 should be small 
or zero. The Higgs sector (2.11c) includes explicit mass terms with real 
parameters Trij^u ^^d' ^^^^ ^ bilinear coupling with complex 
parameter b. Explicit sfermion masses (2. lid) are given by five 3x3 
Hermitian mass-squared matrices mg, m^, m^, m| and nig. The final 
two terms (2. lie and 2.11f) give trilinear couplings between the Higgs and 
squarks or sleptons, with general Yukawa-type complex 3x3 matrices 
Au,Ad,Ae and Cu,Cd,Ce- The latter terms, proportional to the C 
matrices (2.11f), tend to be left out of most low-energy effective models 
because they are suppressed in many SUSY-breaking schemes. 
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>CsuSY-SM is itself straightforward to write down, but lengthy and not 
especially phenomenologically illuminating (see for example Eq. 6.44 of 
Baer & Tata 2006). For the purposes of this general description, the most 
important aspect is that it contains derivatives of the 'superpotential' 



{Yu)ijeabQiH^Uj 

+ 0^d)ijQiHdaDj + {Ye)ijLfHdaEj 



(2.12) 



Here the indices i and j are again generation number, and a and b are 
SU{2)i^ indices. The carets indicate superfields, containing both the SM 
fields (and where applicable, SU{2)i^ doublets) such as ujn^Li, etc. and 
their superpartners ujn^Li^ etc. The terms Dj and Ej are the left 
chiral superfields containing the charge conjugates of the right-handed 
SU{2)i^ singlets: up-type (s)quarks, down-type (s)quarks and (s)electrons, 
respectively. Derivatives of W with respect to its scalar fields give rise to 
all non-gauge interaction terms in /^susy-SM- In this sense, it is some- 
what analogous to the non-gauge part of the scalar potential seen in 
non-supersymmetric quantum field theories. Here it plays a similar role, 
specifying the Higgs potential via the complex parameter /x, as well as 
the Higgs- fermion interactions and ultimately, the fermion masses, via 
the general complex 3x3 Yukawa coupling matrices Y^^, Y^^ and Yg. 

In Eqs. 2.11 and 2.12, only terms which conserve baryon (B) and lep- 
ton number (L) have been retained. B and L are known to be broken 
non-perturbatively in the SM ('t Hooft 1976a, b, see also e.g. the intro- 
ductory remarks of Morrissey et al. 2005). This happens as the B and 
L currents become anomalous (non-conserved) due to a freedom in the 
electroweak vacuum, similar to that discussed for the QCD vacuum in 
the context of the axion in Sect. 2.2. 'Instanton' transitions between 
these vacua can then cause simultaneous violation of B and L by 3 units. 
Nevertheless, B and L are observed to be very good (albeit accidental) 
symmetries of the SM at tree level; the electroweak vacuum doesn't seem 
to have done much tunnelling recently. Exactly why the SM contains 
no B or L-violating operators is an open question. Certain other terms 
which explicitly violate B and L are allowed in Csoft and the superpo- 
tential by all other considerations, but in the interests of constructing a 
model with minimal new interactions, and which agrees with experimental 
constraints, they are generally excluded. 
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Together, >Csusy-SM and the soft terms in Eq. 2.11 give the Lagrangian 
for the the Minimal Supersymmetric Standard Model (MSSM) 



after SUSY breaking but prior to EWSB. The model is minimal in that it 
contains the least additional field content possible for a supersymmetrised 
version of the SM. Note that the only new parameter actually introduced 
in initially supersymmetrising the SM is fi, the Higgs parameter in the 
superpotential. Once we demand a low-energy theory for broken super- 
symmetry, the soft terms introduce more than a hundred additional pa- 
rameters. 

2.4.2 Supersymmetric WIMPs 

In order for a SUSY particle to realistically constitute dark matter, it 
must somehow be stabilised against decay into lighter SM states. The 
most common way this is achieved is to postulate that aside from super- 
symmetry, a discrete Z2 symmetry exists between the SM particles and 
their SUSY partners. The corresponding conserved quantum number is 
known as i?-parity, and has the form 



where s is the particle's spin. All SM particles thus have i?-parity +1 
and all their SUSY partners i?-parity —1. Clearly any theory which si- 
multaneously conserves L and B automatically conserves i?; the MSSM 
superpotential and soft terms (Eqs. 2.11 and 2.12) hence conserve im- 
parity by construction. Note that the converse is emphatically not true: 
i?-parity conservation does not imply B 01 L conservation. For exam- 
ple, the most interesting of the B and L-violating terms discussed in the 
previous subsection (the ones which are arguably valid additions to the 
MSSM Lagrangian), include some which conserve i?-parity. 

Although i?-parity is generally imposed by hand in most SUSY vari- 
ants (either via B and L conservation, or directly), it is not plucked 
entirely from thin air. Some grand unified theories (GUTs) give rise 
to i?-parity conservation naturally in their group and/or representation 
structure; examples are certain supersymmetric 50(10) GUTs (e.g. Lee &; 
Mohapatra 1995). Indeed, i?-parity conservation in a higher theory might 
be part of the reason B and L are approximate symmetries of the SM to- 
day. If i?-parity is conserved, then the lightest SUSY particle (LSP) is 



>Cmssm >Csusy-sm + >Csoft 



(2.13) 



(2.14) 
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absolutely stable. If it is also weakly-interacting and electrically neutral, 
then it is a viable WIMP. 

The quintessential SUSY WIMP is the lightest neutralino, the lightest 
linear combination of the two neutral Higgsinos, the neutral wino and the 
bino, 

XS = V^iHl + V,*iH^, + V,*iW^ + ViiB'', (2.15) 

which mix following EWSB because they share quantum numbers. There 
are thus four neutralinos in the MSSM, corresponding to the different 
linear combinations; only the lightest is ever stable, and it is always this 
particle that is meant when people talk about 'the neutralino' as a WIMP. 
Their masses are given by the eigenvalues of the mass mixing matrix which 
appears in the post-EWSB Lagrangian 



/ n // —9]hL \ 



II 

gvd 

\f2 V2 



gvd 9 vd 
M2 



(2.16) 



in the basis G = {H^, H^,W^ , B^)^ . Here and are the VEVs ob- 
tained by the up- and down-type Higgs during EWSB, and = ^tan^w? 
with the weak mixing angle. The coefficients ^i*/, ^2*/' ^3!' ^1 
individual gaugino/Higgsino components of each mass eigenstate are the 
entries in the inverse of the unitary matrix V which diagonalises TW^o, 



X 



WG, (2.17) 



where X = (Xa^Xb^Xc^Xd)^ • The index / in the V^^, V^^, Vg*^, V^^ refers to 
the column of V which applies to the lightest mass eigenstate. The usual 
labelling for mass eigenstates is to sort X%i Xc^ Xj) increasing mass 
and assign them the numbers from 1 to 4, i.e. Xi? X25X35 Xa- 

Sneutrinos, the spin-0 partners of neutrinos, are also weakly-interacting 
and neutral, so qualify as SUSY WIMPs when the lightest of their number 
is the LSP. Unfortunately, this is not usually the case, because in most 
models one of the sleptons has a slightly lower mass than the lightest sneu- 
trino (Jungman et al. 1996). In any case, constraints on the sneutrino- 
nucleon scattering cross-section from direct detection (Sect. 3.1) all but 
rule it out as the dominant component of dark matter, as the nuclear scat- 
tering cross-section for sneutrinos is even higher than for Dirac neutrinos 
(Falk et al. 1994). 



2.4- Supersymmetry 



33 



It also is worth pointing out that in certain (generahy non-minimal) 
SUSY scenarios, the dark matter particle might not be absolutely stable, 
but metastable on cosmological timescales. In some cases, this removes 
the need for i?-parity conservation. This applies to some types of grav- 
itinos (the SUSY partner of the graviton, the gauge quantum of gravity), 
an example of non-WIMP SUSY dark matter discussed in Sect. 2.5. 

2.4.3 SUSY breaking and parameterisations 

As a low-energy effective theory, the MSSM involves no specification of 
the mechanism by which SUSY is broken. In its full form (Eqs. 2.11- 
2.13), the MSSM possesses a vast number of free parameters that are not 
present in the SM. All but one (/i) come from the soft terms in the SUSY- 
breaking sector (Eq. 2.11). Counting parameters, in the gauge sector we 
have the normal 3 gauge couplings (e,^,^^) and the QCD vacuum angle 
9 (Eq. 2.4) of the SM, as well as 6 gaugino masses (2.11a, 2.11b). One of 
the CP- violating masses can be removed following a field transformation, 
leaving 9 free parameters. From the Higgs sector (2.11c) there are two real 
squared masses Trij^^ and rn^^, and one complex coefficient b. Together 
with the complex parameter and the fact that one degree of freedom 
can again be removed by a field definition, this leaves 5 real parameters. 
In the fermion sector, there are 5 Hermitian 3x3 mass-squared matrices 
(2. lid), as well as 9 general complex 3x3 trilinear coupling matrices 
(2. lie, 2.1 If, 2.12), giving 5 x 9 + 9 x 18 = 207 real parameters. Careful 
treatment of field redefinitions (Baer & Tata 2006) reduces this number 
by 43, leaving 5 + 9 + 207 — 43 = 178 free parameters in the full B and 
L-conserving MSSM, and 178 — 3 x 18 = 124 in the version where the 
'C-terms' (2.11f) are excluded. Given 19 parameters in the SM, there are 
therefore either 159 or 105 more parameters in the full MSSM than the 
SM. 

SUSY-breaking itself is a jungle, and even the zoo of viable breaking 
schemes is vast and varied (see e.g. Luty 2005, for an introduction). Spon- 
taneous SUSY breaking is hard to achieve at tree level, because none of 
the MSSM supermultiplets can easily acquire an appropriate VEV (this is 
related to a particular sum rule for the tree-level masses; see e.g. Martin 
1997; Baer & Tata 2006). The solution seems to be for SUSY to be bro- 
ken in a particle sector which has only very weak or indirect interactions 
with the MSSM, and then communicated to the MSSM by some heavy 
mediator particle X. The mass at which SUSY is broken in the other 
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sector is related to the mass of the mediator by 

^^-msoft-lTeV. (2.18) 

The decoupled matter sector is often referred to by the unfortunate 
sobriquet of 'the hidden sector', which (at least for this author) evokes 
notions of fine-tuning and careful sequestering away of a whole matter 
sector because experimentally, it should be 'seen and not heard'. In reality 
though, most hidden sectors required in SUSY-breaking scenarios have a 
very high characteristic mass scale. We should therefore expect them to 
decouple from the low-energy theory as the influence of the heavy particles 
is integrated out, into mass dimension d > 4 operators suppressed by d — 4 
powers of the large mass scale. In that case, the 'hidden sector' is then 
simply one which is phenomenologically disjoint (effectively hidden) from 
the SM because of its high energy scale. 

Three scenarios have received the bulk of attention to date: grav- 
ity mediation (Chamseddine et al. 1982), gauge mediation (Giudice & 
Rattazzi 1999) and anomaly mediation (Randall & Sundrum 1999). In 
gravity-mediated scenarios, SUSY-breaking occurs in the matter sector 
associated with new physics around the Planck scale Mp. The break- 
ing is hence mediated by particles with Mx ^ Mp, leading to hidden 
sector SUSY-breaking at about 10^^ GeV. The coupling between the hid- 
den and observable sectors is naturally very weak because gravity is so 
weak, and in this case can actually even occur at tree level. This is be- 
cause SUSY-breaking in the gravitational sector actually implies that we 
are talking about supergravity (SUGRA), where SUSY has been made a 
local rather than a global symmetry, and gauged. This happens to mod- 
ify the tree-level mass sum rule in such a way that an appropriate VEV 
is possible to come by from within the MSSM. It is important to note 
that gravity-mediation must always occur at some level if all interactions 
are supersymmetric, the issue is just whether or not it dominates over 
other contributions to SUSY breaking. Aside from string theory, N — S 
SUGRA is also the only known potentially- viable theory of quantum grav- 
ity (Bern et al. 2009). 

Gauge-mediated SUSY breaking (GMSB) takes an entirely orthogonal 
approach, where SUSY breaking is mediated by a particle which has both 
hidden sector and SM gauge couplings. Actually, the mediator gets an 
entire sector of its own in many GMSB scenarios, known as the 'messenger 
sector'; here the hidden sector breaks SUSY, passes it on to the messenger 
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sector, which then communicates it to the observable sector. The media- 
tor conveys SUSY breaking radiatively from the messenger sector to the 
observable sector, coupling to the MSSM only at loop level. In this case 
the mediator can be of almost any mass; if it is light, SUSY breaking in 
the messenger sector can occur at much lower masses than in the gravity- 
mediated scenario. Anomaly-mediated SUSY breaking (AMSB) relies on 
extra dimensions and somewhat contrived brane geometries and dynam- 
ics to suppress gravity mediation to a level where loop contributions from 
a particular anomaly dominate the mediation. 

For the purposes of phenomenology, the full 124 or 178-parameter 
MSSM is not really very practical to work with. Indeed, many of the 
terms should be suppressed anyway in order to agree with experiment: 
the gaugino masses should be small or absent in order to prevent exces- 
sive CP-violation, as should the various CP-violating phases arising from 
complex MSSM parameters. The off-diagonal entries in the mass matri- 
ces are constrained to be small by the level of flavour-changing neutral 
currents (FCNCs) they cause. One approach to SUSY phenomenology 
is thus to work with the effective low-energy theory, but approximate all 
these dangerous terms to zero. One can go a step further and assume 
reality and universality in the mass and trilinear coupling matrices, such 
that each is just a real constant times the identity matrix. One could even 
be so bold as to assume that some of those constants are equal to one 
another. Many different low-energy parameterisations have been used in 
the literature: 6, 7, 8, 19, 24-parameter versions and so on, depending 
upon the individual authors' computational resources and the nuances of 
the specific observable under investigation. 

The alternative is to take one of the scenarios for SUSY breaking, 
which invariably predict universality and reality in the structure of their 
soft terms at the scale of the mediator particles Mx- Investigating the 
phenomenology then requires evolving the soft terms down to the weak 
scale using the renormalisation group equations (RGBs). In the minimal 
supergravity (mSUGRA) model, all scalar masses, gaugino masses and 
trilinear couplings are unified at respective values mo, Mi and Aq, either 

2 

at Mp or more commonly, Mqut- The only remaining non-SM parame- 
ters are tan /3, the ratio of up-type to down-type Higgs VEVs at the weak 
scale, and the sign of the fi parameter. The bilinear coupling B = b/fi 
and the magnitude of /i are set by the condition that SUSY-breaking 
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radiatively generates EWSB. The parameter space is thus^ 

mo, Ml , Aq, tan /3, sgn(//). (2.19) 

The minimal gauge- mediated model (mGMSB) is characterised by the 
parameter set 

A, M, 715, tan/3, sgn(/i), Qrav, (2.20) 

where 

A=^!%^, (2.21) 

\^/ 

with Ms:ygY',mess. the SUSY-breaking scale in the messenger sector and 
(5) the VEV of the scalar part of a gauge singlet superfield S which also 
couples to the messenger sector. M is the mass scale of the messenger 
sector, which consists of generations of vector-like multiplets of quark 
and lepton superfields carrying SM charges. Cgrav is the gravitino mass 
parameter 

^grav — . ' I^Z.ZZ; 

where A%ysY',hidd. is the SUSY-breaking scale in the hidden sector and 
A is the common messenger-sector Yukawa coupling between S and the 
messenger quark and lepton supermultiplets. The trilinear couplings are 
essentially absent in this model, so = = Ag = 0. 
The minimal version of AMSB (mAMSB) is specified by 

mo, ms , tan /3, sgn(/i), (2.23) 

where ms is the gravitino mass, which is also ultimately responsible for 

2 

determining the gaugino masses and the trilinear couplings. The other 
three parameters are defined identically to their mSUGRA counterparts. 



^Strictly speaking, what I describe here is the 'constrained MSSM' (CMSSM). Ac- 
cording to some definitions, the mSUGRA model differs from the CMSSM, in that the 
condition of EWSB does not necessarily have to be imposed, making /i a parameter, B 
is used instead of tan /3, and the relation Aq = B -\- mo is sometimes used to eliminate 
5 as a free parameter (or equivalently, tan/3 in versions where EWSB is imposed). 
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2.4.4 SUSY scanning 

In order to compare any particular version of SUSY with experimental 
data, one generally needs to consider a number of different points in the 
model parameter space, covering a large range of each of the free pa- 
rameters. One then scans over the parameters using some sort of search 
algorithm, evaluating the values of the relevant observables at each point. 
These can then be compared with existing data to determine if the point 
is consistent with experiment, or to make predictions for future experi- 
ments. 

With such high-dimensional parameter spaces, and a very nonlinear 
mapping from some of the model parameters to actual observables like 
cross-sections, it is no surprise that scanning SUSY parameter spaces is 
not a simple exercise. This is especially so if one wants to make meaningful 
statements about which regions of the parameter space are allowed or 
disallowed, and to what degree of statistical significance. Traditional 
analyses simply assigned a status of either 'allowed' or 'excluded' to a 
point if it lay within or beyond a certain confidence level of the measured 
data. These were often based on brute-force random scans, with little 
concern for convergence or any sort of statistical interpretation. The 
modern approach is to consider the full likelihood for each point, and 
analyse the resultant map in terms of either a Bayesian or frequentist 
statistical framework. These scans have often also employed sophisticated 
scanning techniques like Markov Chain Monte Carlos (MCMCs), nested 
sampling or genetic algorithms. The reader is referred to p. 3 of Paper 
V for a detailed background on the differences between the Bayesian and 
frequentist approaches in this context. Details of the differences between 
scanning techniques and a review of significant work in SUSY scanning 
can be found on p. 3 of Paper IV and pp. 4-5 of Paper V. 

In particular, the most well-known products of SUSY scanning are the 
various regions of mSUGRA compatible with the relic density from the 
CMB. These regions, as well as benchmark points drawn from them, have 
been widely used for collider studies and design over the last 20 years. A 
short review of these regions is given on p. 3 of Paper IV. 

2.5 Gravitinos and Axinos 

If supersymmetry is made local and extended to supergravity, the gravi- 
ton acquires a spin-| superpartner, the gravitino. In gravity- mediated 
SUSY-breaking scenarios, the gravitino obtains a similar mass to the 
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other sparticles, i.e. mgoft, whereas in gauge-mediated schemes it has a 
mass of the order of a few keV (Martin 1997). It can therefore be the 
LSP in some cases, the next-to-hghtest SUSY particle (NLSP) in others, 
or just another more massive member of the SUSY spectrum. 

If it is the LSP and i?-parity is conserved, the gravitino is stable, so 
it can constitute a viable but depressing dark matter candidate. This is 
because it interacts only gravitationally with the rest of the spectrum, 
easily fulfilling the collisionless, non-dissipative and electrically-neutral 
requirements of a good dark matter candidate - but the very weakness 
of the interaction means that it is however virtually undetectable. If the 
gravitino is the NLSP, the weakness of its interactions can make it meta- 
stable, potentially leading to non-thermal LSP dark matter production 
by decay. The caveat is that the long lifetime can interfere with later- 
time processes like BBN, depending upon which SM particle the gravitino 
predominantly decays into along with the LSP. The same is also true if the 
gravitino is the LSP, as in this case the weakness of its interactions again 
make the NLSP long-lived. Their long lifetimes mean that gravitinos 
can be viable dark matter candidates even if i?-parity is violated, since 
they would again be meta-stable. Because thermal gravitinos decouple 
from the primordial bath around the Planck mass, their cosmology is 
particularly complicated and uncertain; they could constitute either warm 
or cold dark matter, depending upon their mass and cosmological history. 

If one extends the SM to include the axion (cf. Sect. 2.2) as a solution 
the strong CP problem, then the supersymmetrised version of the theory 
of course also contains a spin-| axino (see e.g. Covi & Kim 2009). Like 
the gravitino, the axino is extremely weakly interacting, so would be 
perfectly viable as dark matter if it were the LSP. In contrast to the 
gravitino though, the axino should at least be detectable indirectly though 
the observation of the axion itself. Because of its ultra weak coupling to 
other particles, if it is the LSP or NLSP the axino might have any of the 
same features and issues as the gravitino with regard to long-lived states. 

2.6 Sterile neutrinos 

The fact that neutrinos have been observed to oscillate indicates that 
they have masses, though experimental limits constrain them to be sub- 
eV. This presents quite a fine-tuning problem, as the masses must be 
many orders of magnitude less than any other known massive particles. 
The favoured mechanism for producing such small masses is known as 
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the see-saw, which operates by introducing heavy singlet (right-handed) 
neutrino states of mass ^M. The observed neutrino masses become in- 
versely proportional to M when the mass matrix is diagonalised. The 
singlet neutrinos do not interact with other particles because they carry 
no charges under any of the SM gauge groups. Being massive fermions, 
the actual neutrino mass eigenstates differ from the weak eigenstates; the 
splitting is sufficiently large that the light neutrinos are almost entirely 
left-handed in character, and the heavy ones almost entirely singlet. This 
makes the heavy states very weakly interacting, or 'sterile', and suffi- 
ciently long-lived to constitute dark matter. 

M can in principle take on almost any value (there are valid natural- 
ness arguments for it being around Mqut, ot for being rather small). If it 
is in the keV mass range, then the sterile neutrino can be a viable warm 
dark matter candidate, with the correct relic abundance (see e.g. Kusenko 
2009, for a recent review). Depending on the production mechanism in 
the early Universe, it is also possible to make sterile neutrinos cold(ish) 
dark matter. 

2.7 Other candidates 

I will not attempt to detail the myriad other viable dark matter candidates 
which exist; this subsection merely provides a summarising 'sound bite' 
on some of the more notable ones, and some references for further reading. 

A range of 'WIMP-derivative' models build upon the basic idea of 
WIMP dark matter, or at least make use of a number of its desirable 
features. WIMPless dark matter (Feng & Kumar 2008; Feng 2010) em- 
ploys the WIMP miracle in a hidden sector, by tuning the WIMP mass 
and interaction strength to achieve the desired relic density, a la Eq. 2.9. 
WIMPzillas (Kolb et al. 1999) are superheavy, weakly-interacting parti- 
cles formed out of thermal equilibrium, arising from e.g. the GMSB mes- 
senger sector. Minimal Dark Matter (Cirelli et al. 2006; Cirelli &; Stru- 
mia 2009) involves the trial-and-error addition of any colourless SU{2)l 
multiplet with a neutral lightest member to the SM, in the search for a 
WIMP. Inelastic Dark Matter (iDM; Smith & Weiner 2001; Tucker-Smith 
& Weiner 2005) and eXciting Dark Matter (XDM; Finkbeiner & Weiner 
2007) are WIMPs with excited states, nowadays often paired with models 
exhibiting Sommerfeld-enhanced annihilation (Arkani-Hamed et al. 2009; 
Nomura & Thaler 2009). These models are all designed to explain specific 
observational anomalies. 
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SIMPs (Starkman et al. 1990) are Strongly Interacting Massive Parti- 
cles which could form colourless bound states (Kang et al. 2008) and hide 
their strong interactions, whilst milli-charged particles (Holdom 1986; 
Davidson et al. 2000) might manage to appear dark because they carry 
only a very small fractional electric charge. Both these options are very 
strongly constrained at the present time (Taoso et al. 2008). 



Chapter 3 



Searches for dark matter 



Dark matter can be sought in a number of complimentary ways. In 
Sects. 3.1-3.5 I give an overview of the main techniques used to search 
for WIMPs, followed in Sect. 3.6 by a quick exposition of other tests, in- 
cluding those for non-WIMP models. It is unlikely that any single search 
technique will ultimately be enough for us to confidently characterise or 
completely exclude a particular dark matter model; combining results 
from different searches is precisely the strategy adopted in Paper IV and 
Paper V, and discussed in Sect. 2.4.4 in the context of SUSY. 

3.1 Direct detection 

Because WIMPs interact with SM particles by the weak force, they should 
have weak-scale scattering cross-sections with normal nuclei. One of the 
most promising ways to detect WIMPs is thus to look for nuclear recoils 
in large- volume target materials on Earth (Goodman & Witten 1985). 

The expected number of WIMP-nucleon scattering events dN per nu- 
clear recoil energy window dE^ is given (Gaitskell 2004) by 



where a is the WIMP-nucleus cross-section, p the local WIMP density, 
the WIMP mass and F the nuclear form factor. The WIMP-nucleus 
reduced mass is /x = (m^mnuc) / (^x ^nuc)- The distribution of WIMPs 
in the halo with velocities v is given by f{v)^ which is integrated over 
all possible velocities. These range from v^i^^E^)^ the minimum velocity 
necessary to produce a recoil of energy £'r, to the halo escape velocity v^^c 




(3.1) 
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Nuclear recoil searches attempt to measure directly. By assuming a 
certain halo model, this then allows a relation to be inferred between 
the WIMP mass and cross-section. The standard halo model assumes 
a Maxwellian velocity distribution with vyims ^ 220kms~^ and a local 
density of 0.3 GeVcm"^; alternative halo models and their influence upon 
direct detection results have also been discussed (e.g. Belh et al. 2002; Stiff 
& Widrow 2003; Green 2007, 2008; Read et al. 2009). 

The cross section has two parts: a spin-independent interaction asi 
between WIMPs and all nucleons, and a spin-dependent component asD 
coupling only to nucleons with net spin. The kinematics of WIMP-nucleon 
collisions mean that both parts are proportional to the square of the 
WIMP-nucleon reduced mass. This dependence is generally suppressed 
in practice however, as direct detection experiments are never sensitive 
to the full WIMP recoil spectrum. The spin-independent part has a fur- 
ther dependence on the square of the nuclear mass. The spin-dependant 
component instead depends upon the nuclear spin J and the spins of the 
individual proton and neutron subsystems in the nucleus (see e.g. Eq. 
2.8 of Paper II, or Jungman et al. 1996; Gondolo et al. 2004; Cerdeho & 
Green 2010). Target nuclei with different isotopic compositions can be 
chosen to optimise an experiment for spin-dependent or spin-independent 
searches. 

Three physical consequences of nuclear recoils are used to search for 
evidence of WIMP scattering. One is ionisation of target atoms caused by 
energy transfer from the recoiling nucleus. Another is fluorescent radia- 
tion given off by electrons of target atoms, as they decay after having been 
excited by energy transfer from the recoiling nucleus. Materials known 
as scintillators, which are transparent at their fluorescent wavelength and 
have a very fast decay time, are used extensively for this purpose; in this 
case, the process is referred to as scintillation. The third technique is to 
measure phonon excitations generated in crystals by the nuclear recoils, 
where minute heat changes corresponding to the absorption of individual 
vibrational quanta are measured. Various combinations of these detection 
methods are employed in different experiments. 

Due to the proper motion of the solar system within the galactic halo, 
there should be a mean net velocity between WIMPs and the Earth. This 
should result in both a diurnal modulation of WIMP-nucleus collisional 
directions, and an annual modulation of the total detection rate due to 
the non-perpendicularity of the solar ecliptic and galactic planes. 

Exclusion limits placed upon the WIMP mass and spin-independent 
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ZEPLIN III (Dec 2008; 
Lebedenko et al. 2009) 

XENON 10 (2007 result, 136 kg d 
exposure; Angle et al. 2008) 

CDMS II (Soudan 2004-2009, Ge 
target; Ahmed et al. 2010) 

SuperCDMS, projected (2 towers, 
Soudan; Akerib et al. 2006) 

SuperCDMS, projected (7 towers, 
SNOLab; Akerib et al. 2006) 

DAMA/LIBRA, with ion 
channelling (5cr; Savage et al. 
2009b, Fig. 11) 
DAMA/LIBRA, no ion 
channelling (5cr; Savage et al. 
2009b, Fig. 10) 

CMSSM Bayesian posterior (68% 
contour; Trotta et al. 2008) 

CMSSM Bayesian posterior (95% 
contour; Trotta et al. 2008) 



Figure 3.1. Exclusion curves, fit regions and SUSY predictions for 
WIMP masses and spin-independent nuclear scattering cross-sections. 
The areas above each curve are excluded by the respective experiments, 
listed to the right. The red regions are claimed detections of an annual 
modulation signal by the DAMA collaboration, under different assump- 
tions about channeling by the detector crystal. The yellow and blue re- 
gions are Bayesian posterior predictions within the CMSSM, taking into 
account all other experimental constraints. 




nuclear scattering cross-sections by various direct detection experiments 
are shown in Fig. 3.1. The plot also shows a theoretical prediction from 
a Bayesian CMSSM parameter scan (blue and yellow solid region; Trotta 
et al. 2008), as well as the area compatible (Savage et al. 2009b, red solid 
region) with the claimed detection of an annual modulation signal by the 
DAMA collaboration (Bernabei et al. 2000, 2008). Other experiments 
appear to disfavour a dark matter interpretation for this signal. However, 
the solid compatible region on this plot has been drawn assuming stan- 
dard WIMP, halo and target models; the annual modulation signal can 
be made marginally compatible with other experiments if these assump- 
tions are relaxed (e.g. Savage et al. 2009a, b; Fairbairn & Schwetz 2009). 
An example of how the compatibility region changes when channelling in 
the detector crystals is included (Savage et al. 2009b) is outlined in red 
dashes. A similar story holds for interpretation of the DAMA signal in 
terms of spin-dependent WIMP scattering (not shown). Predicted sen- 
sitivity curves for upcoming experiments are also shown, showing that a 
substantial part of the CMSSM parameter space should become accessible 
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in the near future. 

3.2 Indirect detection 

If WIMPs are thermal rehcs, they should generically possess weak-scale 
self-annihilation cross-sections. This is true whether the particles are 
Majorana or Dirac. In the Dirac case, being a thermal relic implies that 
there is no net matter-antimatter asymmetry, since any excess of WIMPs 
or anti- WIMPs would cause the relic abundance to be determined by 
the asymmetry rather than thermal production. If dark matter is pro- 
duced non- thermally, it could also have a comparable, or even higher, 
self-annihilation cross-section. Except in the case of a non-thermally 
produced Dirac particle with an initial matter-antimatter asymmetry, 
WIMPs should therefore annihilate at a non- vanishing rate today. 

Indirect detection methods aim to detect the primary or secondary 
products of these annihilations, in the form of photons, neutrinos or other 
cosmic rays. The same techniques can also be used to search for the prod- 
ucts of decaying dark matter, such as the LSP in versions of SUSY which 
weakly violate i?-parity. Searches place constraints on annihilation cross- 
sections in annihilating models (and therefore the relic density, in ones 
where the particle is thermally-produced), and on the lifetime in decay- 
ing models. As a two-body process, the annihilation rate is proportional 
to the square of the dark matter density, whereas the single-body decay 
process is proportional to the first power of the density. 

Whether searching for annihilations or decays, the most promising 
targets are those with large dark matter densities and/or low astrophysical 
backgrounds. The Galactic Centre (GC) would seem the most obvious 
target given its distance and dark matter concentration, but it is also 
one of the most difficult areas to work with because of its complex and 
poorly- understood background (Vitale et al. 2009; Acero et al. 2010), and 
uncertain dark matter proffie (Stoehr et al. 2003; Merritt 2010). Better 
prospects might be had just outside the GC, in a broken annulus which 
excludes the galactic disc (Stoehr et al. 2003; Serpico & Zaharijas 2008). 
Dwarf galaxies are a good option because they are extremely dark-matter- 
dominated, leading to a low background, but have the disadvantage of 
also having rather low predicted fluxes (e.g. Bringmann et al. 2009; Pieri 
et al. 2009; Martinez et al. 2009). Smaller, previously unidentified clumps 
of dark matter with no association to any known astrophysical sources 
(e.g. Green et al. 2004a; Green et al. 2005; Kuhlen et al. 2008; Bringmann 
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2009; Paper III) would make for a spectacular, background-free discovery 
if any were detected, but there is significant theoretical uncertainty in 
their expected number, mass and proximity to Earth. Galaxy clusters 
and the total integrated contribution to the extragalactic background have 
also been considered competitive targets. 

The relative merits of the different targets depend upon whether limits 
are sought on annihilation or decay, and the observation channel. Mas- 
sive charged particles like electrons are attenuated and deflected by mag- 
netic fields, causing them to almost always arrive isotropically. Further 
discussion of the uncertainties in dark matter density profiles and the as- 
trophysical production and propagation of cosmic rays can be found in 
Chapter 4. 

The expected differential primary gamma-ray flux per unit solid angle 
(e.g. Bergstrom et al. 1998) from WIMP annihilations is 



where BF is the boost factor from unresolved substructure in the source, 
/ labels different final states, dN^^ /dE is the differential photon yield 
from any particular final state, cr/ is the cross-section for annihilation 
into that state, v is the WIMP relative velocity, and the integral runs 
over the line of sight to the source. In the absence of any bound states 
(i.e. Sommerfeld enhancements), WIMPs move so slowly that they can 
effectively be considered to collide at rest, allowing ajv to be replaced 
with the velocity- averaged term in the zero- velocity limit, {afv)o. 

In the case of neutralino annihilation, three main channels contribute 
to the gamma-ray spectrum. Through loop processes, annihilation can 
proceed directly into two photons (Bergstrom &; Snellman 1988; Bergstrom 
& Ullio 1997) 
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or into a Z boson and a photon (Ullio & Bergstrom 1998) 



dNl 




dE 



= S{E -m^ + 



(3.4) 



This gives a monochromatic gamma-ray line. Because no other known 
process produces such a line, this would be a smoking gun signal for 
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WIMP dark matter. Unfortunately, the loop suppression means that very 
few models actually have substantial branching fractions into monochro- 
matic photons. A hard spectrum can also be produced by internal brems- 
strahlung (final-state radiation plus virtual internal bremsstrahlung) , gen- 
erated when a photon is emitted from a virtual particle participating in 
the annihilation (Bergstrom 1989; Bringmann et al. 2008). Continuum 
gamma-rays can also be produced by annihilation into quarks, leptons 
and heavy gauge bosons (including the Z from the line), which subse- 
quently decay via tt^ and the emission of bremsstrahlung to softer photons. 

Indirect detection with photons is currently dominated by large air- 
Cerenkov gamma-ray telescopes (ACTs) and the Fermi Large Area Tele- 
scope (LAT; Atwood et al. 2009), a pair-conversion gamma-ray space 
telescope. Limits on the total cross-section from dwarf galaxies (Lom- 
bardi et al. 2009; Essig et al. 2009; Paper IV; Abdo et al. 2010b), the 
isotropic diffuse background (Abdo et al. 2010c; Abazajian et al. 2010) 
and galaxy clusters (Ackermann et al. 2010) are approaching the canon- 
ical thermal cross-section (3 x 10~^^cm^s~^; cf. Sect. 2.3.1), but cannot 
yet conclusively exclude it for any range of WIMP masses. LAT lim- 
its on annihilation into gamma-ray lines also exist (Abdo et al. 2010a). 
Claims have been made from outside the collaboration of excesses in pub- 
lic Fermi data in the inner Galaxy (Goodenough & Hooper 2009; Dobler 
et al. 2010), which some have interpreted as signals of WIMP annihi- 
lation. Their status as true excesses rather than instrumental foibles or 
systematics arising from overly simplified background modelling (e.g. Lin- 
den & Profumo 2010) have yet to be confirmed or disproved by the LAT 
collaboration. 

Gamma-ray and longer-wavelength observations also make it possible 
to hunt for secondary photons from WIMP annihilation or decay. In 
particular, inverse Compton scattering of CMB and interstellar radiation 
field photons by primary leptons injected in annihilations or decays can 
lead to substantial signals in gamma rays and X-rays (Regis & Ullio 2008; 
Cirelh & Strumia 2009; Profumo & Jeltema 2009; Belikov & Hooper 2010; 
Abdo et al. 2010c; Ackermann et al. 2010). Synchrotron emission from 
primary leptons in regions with significant magnetic fields can also be a 
competitive probe (Regis & Ullio 2008; Bertone et al. 2009; Bergstrom 
et al. 2009a), as can modifications of the CMB by particle injection from 
dark matter annihilation at early times (Slatyer et al. 2009; Galli et al. 
2009). 

Recent electron and positron cosmic ray data from Fermi (Abdo et al. 
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2009), the PAMELA satellite (Payload for Anti-Matter Exploration and 
Light-nuclei Astrophysics; Adriani et al. 2009a) the ATIC balloon mission 
(Advanced Thin lonisation Calorimeter; Chang et al. 2008) and the HESS 
ACT (High Energy Stereoscopic System; Aharonian et al. 2008) indicate 
a slight excess in the total number of events at ^10 GeV-10 TeV over that 
expected purely from background. A substantial excess is also seen in the 
fraction of events due to positrons. This has lead to a flood of papers 
explaining the data in terms of dark matter annihilation or decay (e.g. 
Bergstrom et al. 2008; Arkani-Hamed et al. 2009; Nomura & Thaler 2009; 
Chen et al. 2009; Donato et al. 2009). Almost all annihilation scenarios 
require a substantial boost to the annihilation cross-section above that 
expected for a typical thermal relic, either from substructure or specific 
particle models (these are the Sommerfeld-enhanced models mentioned 
in Sect. 2.7). Such boosted models are put under rather severe pressure 
by limits from recent gamma-ray observations (e.g. Abdo et al. 2010b, c; 
Ackermann et al. 2010; Abazajian et al. 2010). The positron excess can 
also be explained in terms of conventional astrophysical sources like pul- 
sars (Profumo 2008; Yiiksel et al. 2009; Hooper et al. 2009) and supernova 
remnants (Blasi &; Serpico 2009; Piran et al. 2009; Fujita et al. 2009), and 
possibly even standard secondary production (Katz et al. 2010). 

In many models, WIMP annihilation produces a substantial num- 
ber of antiprotons. Observations of cosmic ray antiprotons by PAMELA 
(Adriani et al. 2009b) show no excess above the expected astrophysical 
background. This places rather severe limits on the types of models which 
can be invoked to explain the PAMELA positron excess, as they must not 
overproduce antiprotons. Together with the steepness of the rise observed 
in the positron fraction, this suggests that such models should annihilate 
predominantly to lepton-antilepton pairs (so-called 'leptophilic' dark mat- 
ter; Chohs et al. 2009; Bergstrom 2009; Bergstrom et al. 2009b), which is 
somewhat difficult to achieve in the MSSM (Bergstrom et al. 2008). 

Antideuterons are another promising indirect detection channel. Al- 
though predicted yields from dark matter annihilation (Donato et al. 2000; 
Baer & Profumo 2005) and decay (Ibarra & Tran 2009) are certainly 
smaller than those for antiprotons, the very low backgrounds expected 
(Donato et al. 2008) could make antideuterons a useful detection method. 
This is especially true for low WIMP masses, where kinematics would pre- 
vent background antideuteron formation from spallation. Future experi- 
ments, such as the imminent Alpha Magnetic Spectrometer space shuttle 
mission (AMS-02; Choutko & Giovacchini 2008) and the dedicated bal- 
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loon mission GAPS (General AntiParticle Spectrometer; Hailey 2009), 
should make interesting inroads into model parameter spaces. 

In certain models, such as Kaluza-Klein dark matter in Universal Ex- 
tra Dimensions (UED), WIMPs can annihilate directly into neutrinos. 
Other models produce neutrinos in secondary decays and/or cascade in- 
teractions with baryonic matter. Whilst their weak interactions, small 
masses and atmospheric background make neutrinos a difficult indirect 
detection prospect, they have the advantage of pointing directly back 
to a source in the same way photons do. With observations of the GC 
(Bertone et al. 2004) and the inner Galactic halo (Yiiksel et al. 2007), the 
SuperKamiokande neutrino telescope has been used with some success 
to constrain the annihilation cross-section of traditional WIMP models. 
Upcoming cubic-kilometre telescopes such as IceCube should also prove 
powerful enough to essentially rule out Sommerfeld-enhanced leptophilic 
models, through observations of the central regions of our own Galaxy 
(Hisano et al. 2009; Liu et al. 2009; Buckley et al. 2010; Mandal et al. 
2010) or dwarfs (Sandick et al. 2010). 

Indirect detection could also include dark stars and searches for neutri- 
nos from WIMP annihilation in the Sun. Unlike the indirect searches just 
described, stellar probes depend more upon nuclear scattering than the 
annihilation cross-section, so I will describe them separately (in Sects. 3.4 
and 3.5, respectively). 

3.3 Accelerator searches 

Finding dark matter at accelerators is unlikely to be straightforward. The 
most obvious collider WIMP signature is expected to be missing trans- 
verse energy (missing E^)^ which refers to an apparent missing component 
of the total final-state momentum in the direction transverse to a collider 
beam. If the outgoing transverse momenta of a reaction do not sum to 
zero (as they do in the original beam), this can be attributed to the pro- 
duction and escape of a massive particle with a very small interaction 
cross-section with the detector material. The catch is that whilst the 
particle must be stable enough to travel beyond the detector, there is no 
way of to know how stable it is on cosmological timescales. There would 
hence be no way to infer its relic abundance. 

By looking at kinematic endpoints in the momentum distributions of 
observed particles, one can derive the masses of intermediate states in the 
decay chain leading to the missing E't, including the mass of the missing 
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particle itself (Battaglia et al. 2004; White 2007). By carefully examin- 
ing the shapes of the distributions, in principle one can also distinguish 
between different theories giving rise to the same (observable) decay prod- 
ucts. The difficulty with this prospect is that in a hadron collider like the 
LHC (Large Hadron Collider), collisions take place between bound-state 
quarks. It is therefore impossible to ever know exactly where the rest 
frame of each interaction is. Whilst there should be a known distribu- 
tion of rest frames across all collisions, this still introduces a significant 
uncertainty to the distributions for rare processes like those that would 
produce WIMPs, and also makes efforts to disentangle the spins of the 
new particles very difficult (Baltz et al. 2006). 

The main alternative is to attempt to reconstruct the whole underlying 
theory responsible for the new TeV-scale physics, by using other channels 
(than missing E^) to constrain the masses and couplings of other new 
particles. In this way, one might infer the mass, spin and cross-section of 
the dark matter particle(s) without measuring them directly. A more re- 
liable option would be to employ a lepton collider like the proposed ILC 
(International Linear Collider), which should allow direct pair produc- 
tion of WIMPs if they exist in the appropriate mass range. This would 
allow direct measurements of the spin, couplings and mass of the particle, 
though the issue of proving stability would remain (Baltz et al. 2006). 
The couplings could however be measured well enough to determine a 
particle's relic density to a similar level of accuracy as from the CMB 
(Battagha 2009); if the two matched, stability would be strongly imphed. 

Another interesting alternative allowed by a lepton collider would be 
to search for initial state internal bremsstrahlung in dark matter pair 
production (Birkedal et al. 2004). If dark matter were a thermal relic, 
the event rate would be related to the relic density via the produc- 
tion/annihilation cross-section. Some signal should thus be expected for 
any thermally-produced WIMP, assuming a non-vanishing annihilation 
branching fraction into the particular leptons used. 

In general, an unequivocal identification and characterisation of the 
particle responsible for dark matter will require verification from a range 
of different searches. Accelerators will help most in characterising the 
theory to which it belongs, whilst direct and indirect searches should 
provide information about its stability and cosmological abundance. All 
can contribute to determining the mass and couplings. 
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3.4 Dark stars 

If dark matter annihilates at a non-negligible rate today, we might expect 
that the energy from those annihilations could impact nearby baryonic 
systems. This might contribute to heating of the intergalactic medium 
(Ripamonti et al. 2007), reionisation (Natarajan & Schwarz 2008), the 
CMB temperature fluctuations (Slatyer et al. 2009; Galli et al. 2009), or 
the energy budget of stars (Steigman et al. 1978; Bouquet &; Salati 1989; 
Salati & Silk 1989). Stars whose structure or evolution are affected in 
this way are often referred to as 'dark stars'. 

Stars might obtain a substantial amount of dark matter in their cores 
by nuclear scattering, gravitational contraction, or both. In the same 
way that we expect WIMPs to collide with nuclei in direct detection 
apparatuses (Sect. 3.1), we should also expect them to scatter on nuclei 
in astrophysical objects. Thus, in the first case WIMPs could scatter 
on nuclei in stars, lose sufficient energy to become gravitationally bound, 
and return to repeat the process (Press & Spergel 1985; Gould 1987) until 
they either thermalise with the stellar core or annihilate with each other. 
In the second case, steepening of the gravitational potential caused by 
dissipative collapse of a baryonic gas cloud could draw dark matter into a 
star during its formation (Spolyar et al. 2008; Natarajan et al. 2009; Freese 
et al. 2009). In either case, conductive heat transport by WIMP-nucleon 
scattering could in principle also affect the stellar structure or evolution, 
but current direct detection bounds on the scattering cross-sections make 
this unlikely (see e.g. Bottino et al. 2002; Paper II). 

The effects upon different types of stars in a variety of locations have 
been considered in recent years: white dwarfs at the GC (Moskalenko & 
Wai 2007; Hooper et al. 2010) and in globular clusters (Bertone & Fair- 
bairn 2008; McCullough &; Fairbairn 2010), population III stars (Spolyar 
et al. 2008; locco 2008; locco et al. 2008; Natarajan et al. 2009; Freese 
et al. 2008a,b, 2010; Paper VI) and main sequence stars near the GC (Pa- 
per I; Proceeding XIII; Paper II; Casanellas & Lopes 2009). A review of 
the historical and recent developments in this field can be found in Sect. 1 
of Paper II. The two different dark star simulation strategies employed by 
different groups are described in Sect. 2.1 of Paper VI. A technical treat- 
ment of WIMP capture by scattering, and the subsequent impacts upon 
stellar evolution can be found in Sect. 2 of Paper II. For constraining dark 
matter properties, dark stars are most useful for putting limits on either 
nuclear scattering cross-sections or the growth of dark matter halos. 
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3.5 High-energy solar neutrinos 

In the same way that WIMPs would scatter on nuclei in dark stars and 
become gravitationally captured, so they should in the Sun as well. The 
primary difference in this case is that the expected capture rates are far 
lower in the Sun than at the GC or in early proto-halos, mainly due to 
the comparatively low dark matter density in the solar system. Although 
annihilation of WIMPs in the solar core should have essentially no impact 
on the Sun's structure or evolution, some fraction of the annihilation 
energy might escape as neutrinos. 

These (9(GeV) neutrinos would be much more energetic than the 
(9(MeV) solar neutrinos from nuclear fusion (e.g Ahmad et al. 2001), 
so would be clear evidence of WIMP dark matter. The only known solar 
neutrino background at such high energies comes from cosmic ray inter- 
actions with the corona, and is expected to be low (Seckel et al. 1991; 
Ingelman & Thunman 1996). The main background concerns are atmo- 
spheric muons and neutrinos, produced by interactions of cosmic rays 
with the Earth's atmosphere. This background can be mostly avoided 
by triggering only on upwards-going events, at times when a terrestrial 
neutrino telescope is pointing away from the Sun (e.g. Abbasi et al. 2009, 
2010). 

The capture and annihilation of WIMPs in the Sun has been used 
together with telescopes such as SuperKamiokande (Desai et al. 2004), 
AMANDA (Ackermann et al. 2006), ANTARES (Lim et al. 2009) and 
IceCube (Abbasi et al. 2009, 2010; Flacke et al. 2009; Blennow et al. 
2010) to place hmits on WIMP-nucleon scattering cross-sections. Because 
the Sun consists predominantly of hydrogen, the most competitive limits 
from neutrino telescopes are on the spin-dependent cross section. Spin- 
dependent limits from IceCube are already stronger than those from direct 
detection (Abbasi et al. 2009, 2010), and should improve significantly as 
the DeepCore section of the detector is added (Ellis et al. 2010). 

Neutrino telescope limits from the Sun are considerably more model- 
dependent than those from direct detection, as the only way to calculate 
the expected neutrino yields is to have a particular model for the anni- 
hilation branching fractions of the WIMP (see e.g. Wikstrom & Edsjo 
2009). Some channels give rise to more neutrinos than others, producing 
different limits depending upon the assumed annihilation channel. KK 
dark matter for example can annihilate directly into monochromatic neu- 
trinos, whilst neutralinos produce neutrinos from a series of secondary 
interactions between stellar nuclei and daughter quarks and gauge bosons 



52 



Chapter 3. Searches for dark matter 



(e.g. Blennow et al. 2008). The expected signal from neutrino telescopes 
also depends critically on the dark matter velocity distribution in the so- 
lar neighbourhood (e.g. Bruch et al. 2009); capture is far more efficient 
for WIMPs which arrive at the Sun with lower velocities, because they 
need not lose as much energy in the scattering process as faster WIMPs to 
become gravitationally bound. Perturbations of WIMP orbits by planets 
(Gould 1991; Lundberg k Edsjo 2004; Peter 2009) adds further uncer- 
tainty to the expected capture rates for a given cross-section. 

3.6 Searches for other dark matter candidates 

Because axions in electromagnetic fields can convert into photons and 
vice versa (Sect. 2.2), the best places to search for them are in objects 
with strong magnetic fields. The Sun is a good target, as some fraction 
of the photons created in the photosphere will convert to axions as they 
pass through the strong magnetic fields of the chromosphere and corona. 
The CAST (CERN Axion Solar Telescope, Arik et al. 2009) experiment 
searches for these axions by attempting to observe the Sun through an 
opaque screen and a very strong magnet, in the hope that some of the 
solar axions will convert back to photons in the local magnetic field, after 
having passed through the screen (which the solar photons could not). 

A similar idea is pursued in the 'light shining through a wall' exper- 
iments (e.g. Robilliard et al. 2007; Pugnat et al. 2008), where a photon 
source is shone against a barrier after passing through a strong magnetic 
field. A photodetector and a second magnet on the other side of the wall 
look for photons regenerated in the magnetic field from axions that have 
passed through the wall. The same strategy could be used to detect ax- 
ions by looking at background X-ray sources through the Sun (Fairbairn 
et al. 2007). Cavity experiments like ADMX (Axion Dark Matter ex- 
periment, Asztalos et al. 2010) attempt to detect halo axions, by tuning 
the frequency of a magnetic field to the axion mass and searching for mi- 
crowaves from resonant conversion. Astrophysical limits on axions can be 
obtained by considering the maximum amount of energy that they could 
carry out of supernovae and stellar cores without exceeding observed cool- 
ing rates (see e.g. Raffelt 2008, for an up-to-date review). 

One of the more promising ways to search for sterile neutrinos is via 
X-ray line emission produced in the loop decay ^ ^v. Such a signal 
has been claimed in Chandra observations of the ultra-faint dwarf galaxy 
Willman I (Loewenstein & Kusenko 2010), but at less than the 2cr level. 
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Microlensing searches might also stih have some distance to run in de- 
termining the identity of dark matter. This is because ultracompact dark 
matter minihalos (Ricotti & Gould 2009, Paper III) could constitute non- 
baryonic MACHOs. These would apparently have escaped microlensing 
searches because they are slightly more extended than traditional MA- 
CHOs, but might appear in future surveys (Ricotti & Gould 2009). 
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When searching for dark matter, uncertainties in a range of necessary 
input data can influence results. The most relevant of these are the as- 
sumed distribution of dark matter, and the backgrounds to which each 
search is subject. Backgrounds are comparatively well controlled and un- 
derstood in terrestrial direct detection and accelerator experiments, but 
not in indirect searches in astrophysical targets. 

In the types of parameter-scanning exercises discussed in Sect. 2.4.4, 
uncertainties in input data can be fully included in the analysis if they 
have been statistically quantifled. In this case the uncertain quantities 
are referred to as 'nuisance' parameters, because they are parameters one 
is not actually interested in, but that effect results nonetheless. Here I 
give a very brief synopsis of the two most dangerous nuisances in dark 
matter searches. I also provide a brief summary of the most important 
nuisance parameters in particle physics, which necessarily arise from the 
SM. 

4.1 The distribution of dark matter 

Simulations of structure formation using just cold dark matter indicate 
that halos develop self-similar density proflles, approximately following 
an NFW (Navarro, Frenk & White 1996) or Einasto (Navarro et al. 2004) 
proflle. Both these proflles include a steep cusp in the central region. This 
would seem to be at odds with kinematic data, which often indicate that 
the central parts of halos posses smooth cores of dark matter (de Blok 
et al. 2001; Gentile et al. 2004; Del Popolo & Kroupa 2009). New simula- 
tions which include baryons (Governato et al. 2010) show that feedback 
from star formation and supernovae might solve this problem, resulting 
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in cores very much like those observed. This mechanism would not be 
feasible for such low-mass objects as ultra-faint dwarf galaxies, as they 
do not posses enough baryons to form sufficient stars and supernovae to 
make the process efficient. On the other hand, kinematic data from the 
ultra-faint dwarfs is so sparse (e.g. Geha et al. 2009) that they show no 
real preference so far for either cores or cusps anyway, so the smallest 
galaxies may still actually possess the cusps predicted by simulation. 

Because annihilation rates depend on the square of particle densities, 
the predicted rates for indirect detection depend very sensitively upon 
the lower end of the halo mass distribution. This is the source of BF 
in Eq. 3.2, the 'boost factor'. The minimum possible mass of a dark 
matter halo depends upon the nature of the dark matter particle. In 
particular, this is set by the particle's mass and couplings, as well as the 
history of its kinetic decoupling and free streaming in the early Universe 
(Hofmann et al. 2001; Green et al. 2004a, 2005), and how those might 
have been influenced by events like the QCD phase transition (e.g Bring- 
mann 2009). For thermal WIMPs for example, the most likely minimum 
mass is -10"^ (Green et al. 2004a, 2005), though values from 10"^ 
to 10~^ Mq are also possible (Martinez et al. 2009). Smaller halos are 
expected to be far more numerous than larger ones in the picture of hier- 
archical structure formation governed by cold dark matter, although just 
how much more numerous is still a matter of considerable debate between 
A^-body groups (Diemand et al. 2008; Kuhlen et al. 2008; Springel et al. 
2008a,b). 

A standard assumption is that the velocities of WIMPs in halos fol- 
low an isotropic, spherically-symmetric, isothermal (Gaussian) distribu- 
tion, with the width determined by the Keplerian velocity at the solar 
position. In fact none of the details of this assumption turn out to be 
correct if one looks at A^-body simulations: velocity distributions dif- 
fer in the radial and angular directions, vary with galactocentric radius, 
and are neither Gaussian nor spherically-symmetric (Hansen et al. 2006; 
Fairbairn &; Schwetz 2009; Paper II). In any case, the distribution must 
be truncated at the local Galactic escape velocity, which also has some 
uncertainty attached to it. Local inhomogeneities in the phase space dis- 
tribution further complicate matters, due to the presence of streams and 
cool clumps. If this were not enough, recent simulations (Read et al. 2009) 
also suggest that a fraction of Galactic dark matter might exist in a disk 
that co-rotates with the baryonic disk. The presence of such a 'dark disk' 
would substantially boost the bottom end of the velocity distribution, re- 
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Figure 4.1. Propagation volumes for direct observations of cosmic rays. 
Positrons sample only the local environment, as they lose energy very 
efficiently via synchrotron emission, bremsstrahlung and inverse Comp- 
ton scattering, and are easily deflected by magnetic fields because they 
are so light. Antiprotons suffer a similar fate, but are not as strongly af- 
fected because they are much more massive. Gamma rays are neutral and 
massless, so point directly back to their production site. From Bergstrom 
(2009). 

suiting in greater solar capture rates and spectacularly improved limits 
on the nuclear scattering cross-sections (Bruch et al. 2009). Other uncer- 
tainties in the velocity distribution tend to impact direct detection more 
than solar capture, as the former is sensitive to WIMPs of any energy 
above threshold, whereas the Sun only really responds to the low- velocity 
part of the distribution. 

4.2 Production and propagation of cosmic rays 

Cosmic rays (CRs) include protons and heavy nuclei, as well as most 
of the interesting species for indirect detection of dark matter: electrons, 
positrons, antiprotons, antideuterons and gamma rays. CRs are produced 
when particles are accelerated in primary sources such as supernova rem- 
nants, interstellar shocks, pulsars and active galactic nuclei (AGN), and 
in secondary interactions or decay of other CRs. Secondary production 
mechanisms include spallation, radioactive decay of unstable nuclei (in- 
cluding by electron capture), pion production and decay in collisions of 
primary CRs with diffuse gas, inverse Compton scattering of diffuse radi- 
ation, bremsstrahlung and synchrotron emission (Strong et al. 2007). 
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The propagation of CRs depends strongly upon the particle in ques- 
tion, as illustrated in Fig. 4.1. Gamma-rays are neutral and massless, so 
arrive at Earth essentially undeflected and unattenuated. Positrons and 
antiprotons are massive and charged, so suffer deflection in magnetic fields 
and energy losses from bremsstrahlung, synchrotron and inverse Comp- 
ton emission (Moskalenko & Strong 1998; Strong & Moskalenko 1998). At 
the energies relevant to dark matter searches, positrons are much more 
strongly effected than antiprotons. These energy loss processes all pro- 
duce photons, resulting in a diffuse Galactic gamma-ray emission centred 
on the GC and Galactic disk (Strong et al. 2000). Similar processes pro- 
duce gamma rays in external galaxies, especially those actively engaged 
in star formation (Abdo et al. 2010). Unresolved external galaxies and 
AGN could be responsible for the near-isotropic diffuse gamma ray emis- 
sion observed by Fermi (Abdo et al. 2010) and its predecessors. 

Modelling the propagation of CRs essentially requires solving a dif- 
fusion equation in the full phase space of the particles (see e.g. Strong 
et al. 2007), and adding additional terms for particle convection and de- 
cay. The diffusion is sometimes solved for analytically in 1- or 2D (e.g. 
Maurin et al. 2001; Putze et al. 2010), or numerically in 2- or 3D (e.g. 
Strong & Moskalenko 1998; Strong et al. 2007). The numerical approach 
allows the inclusion of much more realistic physics, source distributions 
and geometries. At the present time, the analytic approach is the only 
one which allows a full statistical estimate of the relevant parameters (e.g. 
Putze et al. 2009), and also arguably gives a clearer understanding of the 
physics involved. The parameters of the propagation model can be con- 
strained by comparing with primary-to-secondary ratios such as boron 
to carbon, which are not strongly dependent upon the primary injection 
spectrum. 

Given a set of propagation parameters, a theoretical source distribu- 
tion and injection spectrum, and interstellar gas and magnetic field maps, 
one can then derive a predicted phase space distribution of charged CRs. 
The local form of the distribution can then be used as a background pre- 
diction in e.g. positron or antiproton indirect dark matter searches. To 
obtain a similar estimate of gamma-ray backgrounds, the charged CR 
distribution is combined with the gas map to predict yields from brems- 
strahlung and pion decay, with interstellar radiation maps (and less im- 
portantly, the CMB) to predict inverse Compton yields, and with the 
magnetic field map to produce synchrotron yields. 
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4.3 Standard Model nuisances 

In principle, uncertainties in all 19 SM parameters constitute nuisances 
for SUSY scans, and could impact dark matter analyses. In practice 
though, some parameters dominate over others. The most important 
experimental uncertainties in the SM are those attached to the measured 
masses of the particles. As the heaviest fermions, the masses of the top 
and (to a lesser extent) bottom quarks are the worst constrained of all 
the SM particles except the Higgs. Unfortunately, the heaviest particles 
also have the greatest impacts upon SUSY phenomenology. The masses 
of the Z and W bosons are not substantial nuisances, as they are known 
very accurately (Amsler et al. 2008). 

The strengths of the gauge coupling constants also have a substantial 
effect, and can be difficult to constrain accurately because they run with 
renormalisation group evolution. Only the strong and electromagnetic 
couplings are significant nuisances, as the weak coupling is known as a 
function of the electromagnetic coupling and the masses of the W and Z 
bosons. See e.g. Allanach & Lester (2006), Ruiz de Austri et al. (2006) or 
Trotta et al. (2008) for the importance of considering gauge couplings and 
third generation quark masses as nuisance parameters in SUSY scans. 
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Summary of results 



Paper I, Paper II and Paper VI deal with the theory and detection of 
dark stars, at the GC and in the first dark matter halos. 

In Paper I, we perform a first numerical investigation of the possible 
effects of WIMP dark matter upon main sequence stars. We modify a 
simple static stellar structure code to show that stable stars supported 
by the annihilation of dark matter in their cores have surface temper- 
atures and luminosities that make them resemble protostars (Fig. 5.1). 
The stable solutions we find to the four coupled differential equations of 
stellar structure (see e.g. Kippenhahn &; Weigert 1991; Carrol &; Ostlie 
1996; Stix 2002, or any other standard stellar astrophysics text) track the 
protostellar cooling trajectory known as the Hayashi track. The more 
dark matter is added to the stars, the further their positions are shifted 
in the Hertzsprung-Russell (HR) diagram, back up the Hayashi track and 
away from the standard hydrogen-burning solutions of the zero-age main 
sequence (ZAMS). 

We find that the static code is unable to arrive at stable structures 
across the whole Hayashi track, as it fails to find solutions in the regime 
where neither fusion nor WIMP annihilation dominates the stellar energy 
budget. This can be seen as a series of gaps in the evolutionary tracks 
of Fig. 5.1. We investigate these regions further with a full 'dark' stellar 
evolution code DarkStars (Proceeding XI), and show that the absence 
of solutions in the intermediate regime is simply due to the numerical 
limitations of the static code. 

We go on to investigate the structure and evolution of main sequence 
dark stars in detail in Paper II, using DarkStars. We carefully describe the 
input physics for the DarkStars code, and proceed to catalogue the changes 
in evolutionary paths, core temperatures, core densities and convection 
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Figure 5.1. HR-diagram showing steady- state solutions for main- 
sequence stars powered by differing amounts of dark matter annihila- 
tion. The solid black line indicates the zero-age main sequence (ZAMS), 
where standard hydrogen-burning stars reside. Coloured tracks indicate 
solutions for different stellar masses. Stars situated further along the 
respective tracks (i.e. further from the ZAMS) have been provided with 
more dark matter. The tracks strongly resemble the Hayashi track, along 
which protostars cool as they condense to form main sequence stars. From 
Paper I. 



zones that result from different rates of energy injection by WIMPs, as 
a function of the stehar mass. We show that low-mass stars are gener- 
ally much more strongly affected than high-mass ones, and that WIMP 
annihilation in stellar cores can significantly extend stars' main sequence 
lifetimes. 

We investigate what injection rates could be realised in stars orbiting 
close to the supermassive black hole at the GC, paying attention to stellar 
orbits, dark matter density profiles and velocity distributions. We derive 
new expressions for the capture rate from a truncated Maxwellian distri- 
bution of velocities, and outline an approximate fit to the distribution of 
velocities seen in the Via Lactea (Diemand et al. 2007) simulation. We 
find that low-mass stars on highly elliptical orbits near the GC could ex- 
hibit significant signs of WIMP capture and annihilation (Fig. 5.2). The 
orbits required are very similar to the orbits on which more massive stars 
have already been observed at the GC. The weakest point in this con- 
clusion is the requirement of an adiabatically-contracted NFW density 
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Figure 5.2. WIMP-to- nuclear luminosity ratios achieved by stars on 
orbits with 10-year periods around the Galactic centre. Dark matter 
annihilation can produce up to 100 times as much power as nuclear fusion 
in stars on realistic orbits in our own Galaxy. If the Galactic halo has 
been adiabatically contracted (AC+spike), annihilation can equal nuclear 
fusion in stars on orbits with eccentricities greater than e = 0.9, for masses 
less than about 1.5 M0. If not (NFW+spike), stars of a solar mass or 
less require e > 0.99 to approach break-even between annihilation and 
fusion. The arrow indicates that the IM©, AC+spike curve is expected 
to continue in this direction, but converging stellar models becomes rather 
difficult for such high WIMP luminosities. From Paper II. 



profile for the most marked effects of WIMP annihilation to be realised. 

We show in Paper I, and then even more clearly in Paper II, that dark 
stars resemble protostars because the annihilation of WIMP dark matter 
in their cores causes them to cool and expand. This result is in quali- 
tative agreement with earher analytical estimates (Salati &; Silk 1989), 
though our more accurate numerical calculations show that the stars oc- 
cupy slightly different positions in the HR diagram than simple polytropic 
models would suggest. We conclude that the reason for this behaviour is 
the negative specific heat of a self-gravitating system; when extra energy 
is injected by WIMP annihilation, the self-gravitating body cools and ex- 
pands. Cooling and expansion of a stellar core reduces the rate of nuclear 
fusion. Reduced fusion rates mean that stars' core hydrogen lasts longer. 
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leading to the observed increase in (quasi-)main sequence lifetimes. 

The resemblance of the 'WIMP-burning' tracks to the Hayashi track 
can be understood in terms of two aspects of the physics of WIMP an- 
nihilation in stellar cores. Our appreciation of these two aspects was not 
as well developed at the time of writing of either Paper I or Paper II 
as it has since become, so the following discussion is largely absent from 
those manuscripts. The first aspect is simply that the energy injected by 
WIMPs plays a similar role to the gravitational energy released during 
the contraction of a protostar, as they are both responsible for exactly 
the same type of term in the stellar luminosity equation (one of the four 
equations mentioned above). 

The second and most crucial aspect is that the rate of energy produc- 
tion due to dark matter annihilation is not coupled to the gas pressure or 
temperature in any substantial way; the WIMP distribution, and hence 
annihilation rate, depend only very weakly upon the stellar structure. 
Whereas the fusion rate depends in an essential way on the gas equation 
of state in the stellar core, the dark matter annihilation rate essentially 
depends only on the WIMP density, which is mostly decoupled from ev- 
erything else. Greater energy input from fusion would cause the core to 
cool and expand in a similar manner to what is seen to occur due to 
dark matter annihilation. In this case however, the expansion would have 
the stabilising effect of reducing the fusion yield, and causing the core 
to recontract. On the other hand, structural changes caused by energy 
injection from WIMP annihilations have almost no effect upon the net 
rate of annihilation, so there is no stabilising feedback. This is very much 
akin to the energy released in gravitational contraction; the amount of 
energy released depends only upon the changing gravitational potential, 
not directly upon the thermodynamic gas properties. This is also the 
main reason we see the gaps in tracks in Paper I and Fig. 5.1, and see in 
Paper II that although we can find viable solutions in these regions, such 
stars are only barely stable. 

We turn to dark stars in the early Universe in Paper VI, investigat- 
ing their observability with the upcoming James Webb Space Telescope 
(JWST). We compute model stellar atmospheres and synthetic spectra 
for the dark stars of Spolyar et al. (2009), and feed them through the 
JWST filters to ascertain their detectability as single objects. We find 
that single dark stars in the early Universe will not be detectable by 
JWST except when viewed through a substantial gravitational lens, and 
even then only for certain values of their lifetimes and masses. The ex- 
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Figure 5.3. Left (a): predicted apparent AB magnitudes of dark stars 
at z = 6, as a function the central wavelength observed in JWST broad- 
band filters. Each line corresponds to a different dark star model, where 
models differ in their lifetimes and masses. Models differing in stellar 
mass also necessarily have different intrinsic luminosities, surface tem- 
peratures and surface gravities. Colours refer to surface temperatures 
Teff < 8000 K (red), 8000 K < T^s < 30 000 K (green) and T^s > 30 000 K 
(blue). Dashed lines refer to 5cr and lOcr JWST detection limits after 
3.6 X 10^ s (100 hr) and 10^ s of exposure, respectively. The arrow indicates 
the degree to which curves would be shifted upwards if dark stars were 
observed though the lensing cluster MACS J0717. 5+3745. In this case, 
the coolest dark stars would become individually detectable with JWST. 
Right (b): Colour-colour diagram for dark stars and possible interloper 
populations, based on apparent magnitude differences in different JWST 
filters. The very red spectra of dark stars mean that if they are bright 
enough to be seen at all, dark stars would occupy a rather unique corner 
of the colour-colour plot, and therefore be relatively easy to distinguish 
from other objects. Both figures from Paper VI. 



pected fluxes from different dark stars are shown in comparison to JWST 
detection thresholds in Fig. 5.3a. The scale arrow in this figure indicates 
the amount by which the fluxes would be boosted if viewed through the 
lensing cluster MACS J0717. 5+3745, which provides a magnification of 
/i = 160. Their very red (i.e. cool) spectra mean that if they are visible 
at all, high-redshift dark stars will occupy a rather unique position in the 
colour-colour diagram of Fig. 5.3b, so should be clearly distinguishable 
from interloper objects. 

We also show that if early dark stars live long enough, their longevity 
will combine with their redder spectra to produce a rather peculiar feature 
in the integrated spectra of high-redshift galaxies. The flux contributed 
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by the dark stars will tend to amass in redshift space, producing charac- 
teristic red bumps in the spectra and significantly redder overall galaxy 
colours. If at least 1% of stars in early-type galaxies are indeed dark stars, 
this feature should be detectable in galaxy spectra. 

Akrami et al. (2010); Scott et al. (2010a); Zackrisson et al. (2010); 
Scott & Sivertsson (2009); Fairbairn et al. (2008); Scott et al. (2009b) 

Paper III and Paper IV discuss novel aspects of indirect detection of 
WIMP dark matter using gamma rays, whilst Paper IV and Paper V both 
deal with statistical scanning of SUSY parameter spaces. 

In Paper III, we investigate the prospects for indirect detection from 
a new class of dark matter substructure, proposed by Ricotti &; Gould 
(2009). These ultracompact primordial minihalos would be formed by 
gravitational collapse of large-amplitude, small-scale density fluctuations 
in the early Universe. Such perturbations might have been induced by 
phase transitions such as the QCD confinement transition, or specific late- 
time features in the inflaton potential.^ The resultant minihalos could 
consist of PBHs which have later accreted dark matter, or be made almost 
entirely of dark matter. The first case requires 6 > 30% for the initial PBH 
formation (cf. Sect. 2.1), whereas the second only needs 10~^ ^ ^ ^ 30%, 
making the formation of PBH-free minihalos considerably more likely. 

We compute the gamma-ray fluxes expected from ultracompact mini- 
halos without a central PBH, formed in phase transitions in the early 
Universe. We compare the detection prospects for the Fermi-LAT and ex- 
isting ACTs, showing that they provide largely complementary detection 
capabilities across a range of WIMP masses. We show that single mini- 
halos from the electron-positron annihilation epoch would be eminently 
observable today with existing instruments (Fig. 5.4), and in some cases 
should even have been seen by FermVs immediate predecessor, EGRET 
(the Energetic Gamma Ray Experiment Telescope). We find that miniha- 
los from the QCD confinement phase transition, arguably the best candi- 
date transition for producing sufficient density perturbations to form the 
minihalos, should also be detectable today if any exist within (9(1 pc) of 
Earth. If their population is large enough, it is quite likely that at least one 
such minihalo would be close enough to detect with existing instruments. 
We provide an expression for the cosmological density of ultracompact 

^In the submitted version of Paper III (arXiv:0908.4082vl), we referred to these 
objects as 'PLUMs' (Primordially-Laid Ultracompact Minihalos) so as to distinguish 
them from minihalos produced by the standard spectrum of inflationary perturbations. 
Although that acronym did not ultimately survive to publication, we are pleased to see 
that it has resurfaced in the literature anyway (Lacki & Beacom 2010). 
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Figure 5.4. Integrated fluxes above 100 MeV for primordial ultracom- 
pact minihalos at a distance of d = 100 pc, consisting of WIMPs an- 
nihilating into either bb or ji^ ii~ pairs. Curves are shown for different 
phase transitions giving rise to minihalos, and various degrees of adi- 
abatic contraction. Adiabatically-contracted minihalos are assumed to 
have a fraction F of their mass collapsed into a constant-density baryonic 
core of radius 10~^i?h. Also shown are approximate 5cr, power-law, high- 
latitude, point- source sensitivities for 2 weeks of pointed EGRET and 
one year of all-sky Fermz-LAT observations. Solid limits indicate instru- 
ments' nominal energy ranges. Virtually all types of minihalos produced 
in the e^e~ annihilation epoch would be detectable by Fermi after 1 year. 
Minihalos produced in the QCD confinement phase transition would also 
be detectable if one or more were closer than d ^-j 1 pc. From Paper III. 
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Figure 5.5. CMSSM annihilation cross-sections consistent with all ex- 
perimental constraints, with and without additional constraints from 
Fermi-L AT observations of the dwarf galaxy Segue 1. Scans assume that 
dark matter consists predominantly of neutralinos. Plots show regions 
favoured by existing constraints only (left), and with the addition of 9 
months of Fermi observations towards Segue 1 (middle). Also shown 
is the expected impact upon the favoured masses and cross sections if 
no signal is observed from Segue 1 after 5 years (right). Upper subfig- 
ures give profile likelihoods (yellow and red indicate 68% and 95% confi- 
dence regions respectively) , whereas lower subfigures provide marginalised 
Bayesian posteriors (with 68% and 95% credible regions given by solid 
blue contours). Solid dots give posterior means, and crosses indicate 
best-fit points. The bulk of models disfavoured by observations of Segue 
1 are already strongly disfavoured by other constraints, such as the relic 
density. After 5 years of non-detection, some small inroads will be made 
into interesting parts of the parameter space at very low masses and high 
cross-sections. From Paper IV. 



minihalos, which exhibits a strong dependence upon the spectral index n 
of the density perturbations responsible for their formation. 

In Paper IV, we consider the constraints placed upon the CMSSM by 
Fermi observations of the dwarf galaxy Segue 1. We perform full global 
SUSY fits, and combine them with likelihoods from Fermi data to ascer- 
tain which parts of the CMSSM parameter space are most favoured by 
all current data, and to what degree. Specifically, we include constraints 
from the relic density, electroweak precision observables, the anomalous 
magnetic moment of the muon, rare processes in S-physics and accel- 
erator bounds on sparticle and Higgs masses. We perform full spectral 
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Figure 5.6. Best-fit scalar (mo) and gaugino (Mi) mass parameters 
found by genetic algoritiims (left) and nested sampling (right) in the 
CMSSM. Blue and green shading indicate isolikelihood contours corre- 
sponding to 1- and 2(7 confidence regions in the nested sampling scan. 
Circular, square and triangular markers indicate the best fit in the ge- 
netic scan, the best-fit co-annihilation point in the genetic scan, and the 
best fit in the nested-sampling scan, respectively. The genetic algorithm 
finds more very good fits, but nested sampling maps the regions where 
moderately good fits can be found more completely. From Paper V. 



and spatial fits to the Fermi data, including a full convolution with the 
instrumental energy dispersion and point spread function for all 5.5 x 10^ 
CMSSM models we compute. 

We show that the lack of any gamma-ray signal from Segue 1 dis- 
favours some CMSSM models with very large annihilation cross-sections 
and low neutralino masses. These models all possess cross-sections where 
the neutralino would not be the dominant component of dark matter, so 
are already disfavoured by the relic density constraint (Fig. 5.5). We also 
give a topical discussion of the relative merits of the use of the profile 
likelihood and the Bayesian posterior in SUSY scans. We give examples 
of the additional physical insight that can be gained by considering both 
analysis techniques for any given parameter scan, rather than arguing the 
(rather futile) case that one is somehow more 'correct' than the other. 
Despite the weakness of the constraint drawn on the CMSSM parameter 
space, this paper makes a significant and concrete contribution to the 
literature, as it presents the first inclusion of constraints from indirect 
detection in full statistical SUSY scans. It was also the first dark matter 
paper to use Fermi gamma-ray data, from either within or outside the 
LAT collaboration. 

Paper V explores the scanning algorithms used in statistical analy- 
ses of SUSY parameter spaces. In particular, we test whether genetic 
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algorithms might perform more effective scans in the context of the pro- 
file likelihood than Bayesian techniques, such as nested sampling and 
MCMCs. We implement the genetic algorithm code PIKAIA (Charbon- 
neau 1995) in the SUSY-scanning code SuperBayeS (Trotta et al. 2008), 
using the same experimental data as in Paper IV except for observations 
of Segue 1. 

Our results are somewhat surprising: we find a substantial improve- 
ment in the likelihood of the best-fit point when using genetic algo- 
rithms instead of nested sampling to map the CMSSM likelihood surface 
(Fig. 5.6). We also uncover a section of the stau co-annihilation region 
of the CMSSM at large mo, consistent with all constraints. This sub- 
region appears to have been entirely missed in previous scans. Contrary 
to findings by other authors, often based on Bayesian MCMC methods, 
our best-fit point occurs in the focus-point region, rather than the stau co- 
annihilation region. This particular discrepancy likely has more to do with 
differences in physics codes than the scanning algorithm itself. Nonethe- 
less, our direct comparison with nested sampling in SuperBayeS indicates 
that all groups using Bayesian techniques to perform SUSY scans in the 
context of the profile likelihood should have significant concerns about 
whether the convergence of their scans is sufficient for a non-Bayesian 
analysis. 



5.1 Outlook 

This thesis presents a number of novel approaches and techniques for 
discovering or constraining the nature of dark matter. Its results also 
provide some initial applications of those techniques to constraining dark 
matter models, especially SUSY in the form of the CMSSM. 

A myriad of extensions to the studies described here could be imag- 
ined, from hybrid scanning algorithms, additional observables and nui- 
sances in SUSY scans, to post-main sequence evolution of dark stars and 
more detailed modelling of the halos in which they reside. The scan- 
ning techniques of Paper IV and Paper V, along with the improved solar 
abundances of Paper VII, Paper VIII and Paper IX, should help sharpen 
our understanding of the limits on nuclear scattering cross sections pro- 
vided by neutrino telescopes. Together with an analysis of unidentified 
point sources in Fermi data, the cosmological density of ultracompact 
minihalos and flux predictions derived in Paper III would also be useful 
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for constraining the spectrum of perturbations produced in various phase 
transitions in the early Universe. 

Together, Paper I, Paper II and Paper VI indicate that the prospects 
for detecting dark stars or using them to constrain nuclear scattering 
cross-sections are not fantastic, although results are promising enough not 
to write off the possibility. Dark stars in the early universe are probably 
more likely to exist and show effects of WIMP annihilation, but those 
in our own Galaxy might be easier to detect and investigate in detail. 
Based on the results of Paper IV and other careful analyses of Fermi data 
(Abdo et al. 2010b, c; Ackermann et al. 2010), it would seem that indirect 
detection is also unlikely to provide a reliable discovery or exclusion of 
standard WIMP dark matter in the near future. This statement would 
not apply if we happen to be very lucky with respect to the identity 
or early history of dark matter; if it were of a Sommerfeld-enhanced or 
non-thermal variety, or the QCD phase transition in the early Universe 
substantially increased the substructure boost factor through the creation 
of ultracompact minihalos, the picture for indirect would be rather more 
optimistic. 

Notwithstanding these rather fine-tuned scenarios, it would seem that 
the best hope for dark matter detection in the next few years is probably 
at the LHC, and in direct detection experiments. In this case, analyses 
like those of Paper IV and Paper V will become especially important 
as we try to extract information about dark matter from LHC data, and 
cross-correlate it with dedicated dark matter probes. A credible discovery 
of dark matter is unlikely to to be achieved through a single detection 
channel, much less a convincing identification in terms of the underlying 
quantum field theory. 
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